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1 

Dynamical mean-field theories of 
correlation and disorder 



1.1 Mott transitions in clean and disordered systems 



It is this fascination with the local and with the failures, not successes, of band theory, 
which... contradicted the assumptions of the time... 



P. W. Anderson, Nobel Prize Lecture, 1979 



1.1.1 Introduction 



The Dynamical Mean Field Theory (DMFT) of interacting lattice fermions is perhaps 
best described as the optimal description of these systems which takes into account 
only local (on-site) correlation effects (Pruschke et ai, 1995 Georges et al., 1996). 



Technically, this is implemented through the assumption of a single-electron self-energy 
which, in the lattice translationally invariant case, is independent of the site in real 
space or of wave vector in reciprocal space and therefore is a function of frequency 
only 

^ or i;(fc,w) ^ 

This kind of approximation surely leaves out inter-site correlations. However, local ap- 
proaches to strong correlations have a long history and several semi-phenomenological 
descriptions of classes of compounds within this framework have been shown to be con- 
sistent with their physical properties. A particularly well-studied example are heavy 



fermion systems (Stewart, 1984[ Grewe and Steglich, 1991). Many properties of heavy 



fermion materials suggest that local correlations are sufficient for a good description, 
the large carrier effective mass (as derived from the magnetic susceptibility or the 
Sommerfeld specific heat linear coefficient) being the best known but by no means the 
only ones. Transport properties such as DC and optical conductivity and ultrasound 



attenuation can also be understood with the same assumptions (Varma, 1985). Other 
examples of systems well described by a local approach include systems close to a Mott 
metal-insulator transition, in particular, the vicinity of the finite-temperature critical 



end-point ( 


Rozenberg et al., 1999 


Kothar et al., 2000 




Kothar et al, 2002 


Limelette 


et al, 2003a 


Limelette et al., 2003 & Kagawa et al, 200S 


1 . A particularly striking con- 



sequence of a local self-energy is the cancellation of many-body renormalizations in 
the ratio of the coefficient of the T^-term of the resistivity to the square of the specific 



heat coefficient, usually called the Kadowaki- Woods ratio (Kadowaki and Woods, 1986 



Miyake et al., 1989). Recently, it has been shown that, when materials-specific effects 
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(such as carrier density, density of states and Fermi velocity values) are properly taken 
into account, then the Kadowaki- Woods ratio appears to be universal across a much 
wider range of compounds, including, besides heavy fermion systems, organic charge- 



transfer salts, transition- metal oxides, and transition metals ( Jacko et al., 20091. Thus, 
a local approach to strong correlations seems to be much more generally valid than 
initially thought. 

Several starting points lead to theories that ultimately predict a self-energy of this 



form, most notably descriptions based on the Gutzwiller wave-function (Gutzwiller, 



1963 


Gutzwiller, 1965 


Vollhardt, 1984 


and Lee, 1987 


). However, these theories 



beyond a local self-energy, e. g., inelastic scattering effects are not included, higher- 
energy incoherent features are absent, among others. A description which incorporates 
all possible local effects in a fully self-consistent fashion is provided by DMFT. 

Historically, DMFT was proposed by a recourse to the infinite-dimensional limit 



of lattice systems (Metzner and Vollhardt, 1989). Indeed, when appropriate rescaling 



of parameters is done (as is usual when considering this limit), the theory remains 
meaningful and non-trivial as c? — >■ cx) and the self-energy becomes completely local. 
The reader can find many alternative derivations of DMFT in this limit in the review 



(Georges et al., 1996). Alternatively, DMFT can also be viewed as the best local 
description of three-dimensional systems. The focus here will not be to derive the 
theory by resorting to the infinite-dimensional limit but rather to highlight the physical 
content of a local description of correlation effects. This is specially important since we 
will later explore other, more general local theories which do not become exact in any 
particular limit but which inherit the insights gained from DMFT. We will therefore 
focus mostly on a Bethe lattice, which is most transparent and lends itself particularly 
well to generalizations to the disordered case. Furthermore, we will highlight the key 
physical assumptions involved, which are kept in the other approximations. 



1.1.2 The clean case 

Consider for concreteness the Hubbard model with only nearest-neighbor hopping on 
a lattice with finite coordination z, in usual notation. 



H 



(1.1) 



We focus on a particular site, call it j, which in the clean case can be any site. The 
effective dynamics of this site alone can be obtained by integrating out all the other 
sites. This is no longer a Hamiltonian dynamics and the procedure requires an action 
description, which we will write in imaginary time 
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r0 

Seff (j) = / drc]^ (t) {dr - Ai) Cja (t) 



+ U / drnj^ (r) n,; (r) . (1.2) 
Jo 

The second term above comes from integrating out the other sites, the first and third 
ones being the local contributions, already present before the integration. The thing 
to note here is the fact that, in general, the integration over interacting sites gener- 
ates other higher-order terms, involving four and more fermionic fields. In the high- 
dimensional limit or in DMFT in general, these higher-order terms are absent or ne- 
glected. It is clear that this means that only single-particle inter-site correlations are 
kept in this limit/approximation and this is precisely what is encoded in the second, 
retarded term in Eq. ( |1.2[ ) . The "hybridization function" A (r) describes the "leaking" 
of electrons in and out of site j. It can be written as 

A(r)=i2 ^ a'iilir), (1.3) 

/,m— 1 

where (r) is the Green's function for propagation from site m to site I in a lattice 
from which site j has been removed (hence the superscript (j)) 

G[^ir) = -{T[c^Ar)clAO)]f\ (1-4) 

and the sums extend over the z nearest-neighbors of site j. Let us not dwell on how 
this is calculated for now. 

The action (1.2) is equivalent to the one of an Anderson single-impurity problem 



(Anderson, 1961), whose Hamiltonian is 



(^«Lcj- + H.c.^ +Uc]^c,^c]^c,^, (1.5) 



provided we choose Ek and Vk above in such a way that the Fourier transform, in 
Matsubara frequency space, of the hybridization function in Eq. (1.3 1 is such that 

^(-") = 1E^^^- (1-6) 

Ns ^ iw„ - Ek 

k 

This equivalence proves to be extremely useful since the well-studied behavior of the 



Anderson single-impurity problem serves as a guide to physical insight (Georges and 
Kotliar, 1992) . 
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Suppose now that we can somehow find the full interacting Green's function of the 



system described by the action of Eq. ( 1.2 1 



G,, (r) = -(T [c,. {t)cI (0)]) 



(1.7) 



where the subscript ef f emphasizes that it is to be calculated under the dynamics 



dictated by (1.2). We can repackage our ignorance about this function by defining a 



self-energy S (iw„) such that 



1 



+ fl- A [iLOn) - S {iUn) 



lUJ, 



(1.8) 



Having quantified the local dynamics, we now need to bring in information from the 
rest of the lattice. In principle, this is quite straightforward: since only local correlations 
are included, S (iaj„) is also the self- energy for generic lattice propagation 



G (fe, iujn) = 



1 



iUn - efc + /i - S (iUJn) ' 



(1.9) 



where 6^ is the non-interacting dispersion. From this expression, we can obtain the 



Green's function with one site removed from Eq. ( 1.4 1 and from that the hybridization 



function (1.3), which closes the self-consistency loop. This procedure to get A (ja;„) 
from G{k,iujn) is described, e. g., in the review (Georges et ai, 1996). However, 
we will proceed in the simpler and more illuminating case of the Bethe lattice with 
coordination z, for which 



(1.10) 



since the removal of site j completely disconnects the two branches that start at the 
nearest-neighbors I and m, if ^ 7^ to. Thus, 



Mr)^t^Y.G\r{r) 

1=1 

- .t^G« (r) , 



(1.11) 



(1.12) 



since all sites are equivalent. We can now take the limit z — )• 00, noting that, in this 
limit, the removal of one nearest-neighbor site (j) is irrelevant for the local propagation 
at I and that an appropriate rescaling, namely zt^ t^ , is necessary 



A (r) = PGu (r) = PG,, (t) 



(1.13) 



since all sites are equivalent. This is the self-consistency condition in this case: the 
solution of the problem is the one for which, if we plug in the local Green's function 



Gjj (t) in Eq. ( 1.13), then insert it into the action ( 1.2 1 and find the expectation value 



in Eq. (1.7), we get back Gjj (r). 
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A more physical alternative route and one which is not restricted to the Bethe 



lattice is to note that the local Green's function obtained from (1.91 must coincide with 



the one defined in (1.8 I 



1 



* k 

Po (e) 



iuJn + /i — A (iuJn) — S {iOJn) 



= / de 



iUJn — £ + /i — S (iuJn) 



(1.14) 



Here, po (s) is the bare density of states generated from . The two procedures can 
be shown to be equivalent and establish the necessary self-consistency condition. 
It is possible to extend the above analysis to two-particle correlation functions 



(Georges et al, 19961 but we will not delve into it. It is, however, worthwhile to notice 



that the current-current correlation function, which provides the conductivity through 
the Kubo formula acquires no vertex correction within DMFT and is given by a simple 
bubble of renormalized single-particle Green's functions. 

It is perhaps wise to highlight yet again what the key assumptions of this approach 

are: 

1 . The effects of interactions that are included are on-site only, or equivalently, the 
self-energy is purely local. 

2. Different sites "know" about each other through single-particle processes only, see 
Fig-inja). 

3. The local dynamics, as dictated by the local effective action and usually encoded 
in the local Green's function, must coincide with the local dynamics as derived 
from the lattice propagation. 

The DMFT approach, like the original Bragg- Williams (BW) mean-field theory of 
magnetism (Goldenfeld, 1992), focuses on a single lattice site, but replaces its envi- 



ronment by a self-consistently determined "effective medium" (Georges et al., 1996) 



Unlike the BW theory, the effect of the environment cannot be captured by a static 
external field, but must be encoded in a full complex function A (ia;„), which contains 
information about the dynamics of an electron moving in and out of the given site. 
The calculation then reduces to solving an appropriate quantum impurity problem. 



Eq. (1.2), supplemented by an additional self-consistency condition, Eq. (1.14), that 



ultimately determines this hybridization function A(iw„). 

The approach has been very successful in examining the vicinity of the Mott tran- 
sition in clean systems, in which it has met spectacular success in elucidating various 



properties of several transition metal oxides (Georges et al, 1996), heavy fermion 



systems, and even Kondo insulators (Rozenberg et al., 19961 



1.1.3 The clean Mott transition 

The Mott transition in a single-band Hubbard model can be regarded as a prototype 
for a interaction-driven metal-insulator transition, a phenomenon with plausible rel- 
evance to many physical systems of current interest. Its basic mechanism has been 



correctly understood for more than fifty years (Mott, 1949), yet the precise nature of 



this phase transition has long remained controversial and ill-understood. Part of the 
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Fig. 1.1 Schematic view of all DMFT-inspired theories of disordered strongly correlated sys- 
tems, (a) All these theories share the same ingredient: the local site "knows" about the other 



sites through the self-consistently determined hybridization function A{uj) (see Eq. (1.2|), 



which acts as an order parameter. They differ in the level of description of this order param- 
eter, which is in general a random quantity (panels b and c). Both DMFT and the Typical 
Medium Theory (TMT) replace the actual local realization of this random variable by a fixed, 
non-random function, namely, the average and the typical values of A (cj), respectively. These 
two functions are very similar at weak disorder (panel b), but become increasingly different 
as the disorder grows (panel c). The statistical dynamical mean field theory (statDMFT), by 
contrast, retains the full spectrum of spatial fiuctuations, since each site "sees" the actual 
local realization (lu). 

confusion stems from the fact that at low temperatures the Mott insulator is typically 



unstable to antiferromagnetic ordering, leading many authors (Slater, 1951) to focus 



on magnetism as a proposed driving force. The shortcoming of this view was most lu- 



cidly emphasized by Anderson ( Anderson, 1978 ) who stressed that the Mott insulating 



state persists well above the Nel temperature. It is thus transmutation of conduction 
electrons into local magnetic moments - not the long range magnetic ordering - that 



Mott transitions in clean and disordered systems 7 



should be regarded as the fundamental physical process behind the Mott transition. 
The two phenomena can be most clearly separated in systems where the tendency 
for magnetic ordering can be appreciably weakened due to frustration effects, such 
as often found in orbitally degenerate transition metal oxides. Here, the competition 
between antiferromagnetic superexchange and ferromagnetic tendencies due to Hund's 
rule couplings typically lead to large cancellations, resulting in very weak magnetic 
correlations in the paramagnetic phase. 

From the theoretical point of view, this situation can be most clearly formulated by 
focusing on the "maximally frustrated" Hubbard model, with infinite-range hopping of 



random sign (Georges et al, 1996). In this model the magnetic frustration is so strong 
as to completely suppress any magnetic ordering, while the DMFT approximation 
becomes exact, allowing precise and detailed characterization of such an "ideal" Mott 
transition within the paramagnetic phase. In the following we briefly describe the main 
features of the resulting DMFT picture of the bandwidth-driven Mott transition. 

Critical behavior and mass divergence at T = 0. Within DMFT, the critical regime 



between the Fermi liquid metal and the Mott insulator features (Georges et al, 1996) 
a finite temperature coexistence region and a first-order transition line ending at the 
critical end-point at T = Tc. At T — 0, however, the metallic solution is the stable 
(lower energy) one throughout the coexistence regime. It is characterized by heavy 
quasiparticles (QP) with an effective mass that diverges as the transition is approached. 

In general, the effective mass is evaluated from the single-particle self-energy using 
the expression 



l-^S'(k,a;) 



^ l + fMS'(k,..) 



(1.15) 

k — kp ,0;— 



fe dk' 

where kp denotes the Fermi momentum, and S' is the real part of the self-energy 
S(k, oj). The QP weight, on the other hand, is defined by 



(1.16) 

k—kp ,0;— 



Within DMFT, E(k,w) = I](w) is momentum independent, and m* /m = Z ^. Note 
that, since generally -^T,' {uj)\^_^ < 0, the interactions increase the effective mass. The 
actual divergence is obtained only if the quantity A = — g^S'(a;)|^_p itself diverges. 
This scenario is realized, for example, in the Brinkmann-Rice theory of the Mott 
transition, as well as in the more recent DMFT solution. Since the QP weight is 
simply = m*/m, it must diverge at the same place as m* does. 

We should emphasize that this result is exact within the DMFT approach and is 
an excellent approximation for many Mott compounds where magnetic frustration is 
sufficiently strong. To put this result in perspective, we contrast it with a popular 
but uncontrolled weak-coupling approach, based on the so-called "on-shell approxima- 



tion" (Ting et al., 1975) for the effective mass of the correlated electron gas. Here, an 



approximate expression for the effective mass is proposed 
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where is the unrenormahzed band dispersion. When this approximation is apphed to 
the low-density electron gas within the Random Phase Approximation (RPA) scheme 



( Zhang and Das Sarma, 2005 ) , one finds that the efi^ective mass diverges before the QP 
weight Z vanishes. This result seems quite pathological, since the natural interpreta- 
tion of the efi^ective mass divergence is the localization of itinerant electrons, where 
one also expects the breakdown of the quasiparticle picture. 

To benchmark the validity of the proposed "on-shell approximation", we apply it 
to the maximally frustrated Hubbard model, where DMFT provides us with an exact 
result for both the effective mass and the full self-energy. In this case E'(k, cj = i^k) = 



T,'{uj=^u.). Noting that 

m d 



k—ki 



0=0 



we get the "on-shell" result 



m* _ 1 



0=0 



1 

1-A' 



As we can see, this expression is equivalent to the exact expression m* /m = 1 + A, 
only to leading order, i.e. for (1 — <C 1. On the other hand, the positive quantity 
A is expected to grow with the interaction. As long as it is finite, neither will the 
properly defined effective mass m*/m, nor the inverse QP weight ever diverge. 
In contrast, if one uses the "on-shell" expression, then the effective mass will blow 
up as soon as ^4 = 1, and this will happen at some point in any approximation 
where A grows with the interaction. However, as we can see, this will not lead to the 
divergence of the inverse QP weight Z~^. What we can see from these expressions is 
that the essence of the "on-shell" approximation is simply to linearize the expression for 
(m*/™)^^, by expanding it in the quantity A = —^'S'{uj)\^_^ . Instead of appearing 
in the numerator of the effective mass expression, it now enters the denominator, 
leading to an unphysical effective mass divergence. This example provides a perfect 
illustration of how dangerous it is to indiscriminately apply weak-coupling results to 
non-perturbative phenomena near the Mott transition. It also shows how DMFT not 
only correctly captures the essence of strong correlations, but also provides a simple 
and transparent insight into their physical content. 

Quantum-critical behavior at T > Tc- Many systems close to the metal-insulator 
transition (MIT) often display surprisingly similar transport features in the high tem- 
perature regime Here, the family of resistivity curves typically assumes a characteristic 
"fan-shaped" form, reflecting a gradual crossover from metallic to insulating transport. 
At the highest temperatures the resistivity depends only weakly on the control param- 
eter (concentration of charge carriers or pressure) while as T is lowered, the system 
seems to "make up its mind" and rapidly converges towards either a metallic or an in- 
sulating state. Since temperature acts as a natural cutoff scale for the metal-insulator 
transition, such behavior is precisely what one expects for quantum criticality. In some 



cases (Abrahams et al., 2001), the entire family of curves displays beautiful scaling 
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Fig. 1.2 (a) DMFT resistivity curves as functions of tiie temperature along different trajec- 



tories across the Mott transition in a half-filled Hubbard model (Terletska et al, 20111. (b) 



Resistivity scaling displaying remarkable "mirror symmetry" (Dobrosavljevic et al, 19971 



behavior, with a remarkable "mirror symmetry" of the relevant scaling functions (Do- 



brosavljevic et al., 1997). But under which microscopic conditions should one expect 
such scaling phenomenology? Should one expect similar or very different transport 
phenomenology in the Mott picture? Is the paradigm of quantum criticality even a 
useful language to describe high temperature transport around the Mott point? 

Somewhat surprisingly, most DMFT studies of the Mott transition focused on the 
lowest temperature regime, paying little attention to the high temperature crossover 
regime relevant to many experiments. On the other hand, it is well known that at 
very low temperatures T < Tc ^ O.OSTp, this model features a first order metal- 
insulator transition terminating at the critical end-point Tc (Fig. 2), very similar to 
the familiar liquid-gas transition. For T > T^^, however, different crossover regimes have 



been tentatively identified (Georges et al., 1996) but they have not been studied in any 



appreciable detail. The fact that the first order coexistence region is restricted to such 
very low temperatures provides strong motivation to examine the high temperature 
crossover region from the perspective of "hidden quantum criticality" . In other words, 
it is very plausible that the presence of a coexistence dome at T < <C Tp, an 
effect with very small energy scale, is not likely to influence the behavior at much 
higher temperatures T ^ Tc. In this high temperature regime smooth crossover is 
found, which may display behavior consistent with the presence of a "hidden" quantum 
critical (QC) point at T = 0. To test this idea, very recent work (Terletska et al., 2011 ) 



utilized standard scaling methods appropriate for quantum criticality and computed 
the resistivity curves along judiciously chosen trajectories respecting the symmetries of 
the problem. Characteristic scaling behavior for the entire family of resistivity curves 
has been identified, with the corresponding beta-function displaying striking "mirror 
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symmetry" consistent with experiments. These findings provide compelling arguments 
in support of the suggestion that finite temperature behavior in many Mott systems 



should be interpreted from the perspective of quantum criticality ( Panagopoulos and 



Dobrosavljevic, 2005 ). It should stressed, however, that the DMFT solution in question 
does not contain any physical processes associated with approach to magnetic or charge 
ordering. If quantum criticality is indeed at play here, it has a fundamentally different 
nature, one that is associated with the destruction of a Fermi liquid without the aid 
of any static symmetry breaking pattern - in dramatic contrast to most known critical 
phenomena. 

1.1.4 The disordered case 

Let us now proceed to write down the equations for the dynamical mean field theory 



description of a disordered system ( 


Janis and Vollhardt, 1992 


Janis et al, 1993 


Do- 


brosavljevic and Kotliar, 1993 Dobrosavljevic and Kotliar, 1994 


). We will focus on 



the case of diagonal site disorder, which for the Hubbard model reads 



H = - 



E 



(1.18) 



where Sj are assumed to be independent random variables drawn from a given distri- 
bution P (e) whose strength is W (say, a uniform distribution from —W/2 to W/2). 
Focusing once more on the local dynamics, it is clear it must now be dictated by 

Seff U) = E / (^) i^r + £j - Cja (t) 

+ E / / dr'4(r)A(r-r')c,.(r') 
„ Jo Jo 



+ U / drrij^ (r) n,^ (t) . 
Jo 



(1.19) 



Notice that the effective action is now different for different sites because of the £j term. 
However, the hybridization function A (t) is not site-dependent. This is motivated 
again by the infinite-dimensional limit. Indeed, if the number of nearest-neighbors is 



infinite, the sum in Eq. (1.3) is effectively an averaging procedure over all possible 
realizations of the Green's function. In fact, in the infinite-dimensional Bethe lattice, 
Eq. becomes 



A (t) = zt^ 



z 

L 1=1 



-^t'Gu (r). 



(1.20) 
(1.21) 



where the overbar denotes average over quenched disorder. Again, the effect of remov- 
ing one nearest-neighbor is negligible in this case. Thus, the hybridization function is 
proportional to the average local Green's function (see Fig. |1.1[ 3 and c). Solving 



Mott transitions in clean and disordered systems 11 



the DMFT equations in the disordered case then entails solving an ensemble of single- 
impurity problems as in Eq. (1.191, one for each value of Sj, and finding for each of 



them the local Green's function and self-energy (which are now also site-dependent) 



Gnir) = -{T [c,. (r) c], (0)] ) 



IUJ„ 



(1.22) 
(1.23) 



The self-consistency can then be written, for the Bethe lattice, as (cf. Eqs. ( 1.11|1.13 1) 

P(e) 



A (iw„) = t Gjj (iun) = t 



de 



iujn - e + M - A (iujn) - S [e, iuJn] ' 



(1.24) 



where we have slightly modified the notation in order to show that the denominator 
on the right-hand side depends on the site-energy e both explicitly and implicitly 
through the self-energy S [e, ia;„]. It is obvious that the above procedure reduces to 



the original DMFT in the clean case (cf. Eq. (1.13)) but what does it reduce to in the 
non-interacting, disordered case? 

It turns out that the treatment of disordered non-interacting systems obtained 
from these equations is equivalent to the so-called Coherent Potential Approximations 



(CPA) (Elliott et al., 1974 Economou, 2006), which is known to become exact in 
infinite dimensions (Vlaming and VoUhardt, 1992 1. The CPA equations are usually 



obtained through a strategy that consists in replacing the effects of scattering off 
the exact disorder potential by an effective average medium. Formally, one writes the 
average Green's function in terms of an average medium self-energy Sam {i^n) 
frequency-dependent complex quantity) 



G{k,k',iu}n) 



->k.k' 



(1.25) 



where again is the clean non-interacting dispersion. S^Af (i^n) is calculated by 
replacing the average medium by the exact potential (as defined by the actual values 
of Ej) at a single generic site (while keeping it at the other sites) and imposing that 
the difference between the exact and the average scattering t-matrices vanishes on the 



average (Elliott et al., 1974 Economou, 20061. A similar, effective medium approach. 



incidentally, can be used to derive the DMFT of clean interacting systems (Georges 
et al., 19961, so it is no surprise that one recovers CPA in this case. In the generic 
case of disordered interacting systems, Sam (i'-^n) is obtained from the local part of 
the average Green's function 



Gjj (iUJn) 



Pie) 



de 

lujn - e + fi- A [lujn) - S [e, iw„J 

±y I . 

Ns ^ lUJn - £fc + A* - ^AM («W„) 



(1.26) 



One may wonder what is the form of the DMFT self-consistency for generic disor- 
dered interacting systems, beyond the Bethe lattice case, in other words, the analogue 
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of Eq. (1.14). This is most easily done through the analogy with CPA. Once we have 
the local self-energy for every value of Sj, T,j {iUn), for a given A Eqs. ( 1.22|1.23 ), 

we first find the average medium self-energy T,am («w„) through Eq. (1.261. We then 



note that the average local Green's function within CPA is also given by 



Gjj (iUJn) 



1 



iuJn + M - A (iuJn) - [i^n) ' 



(1.27) 



since this is what you get if you replace the actual scattering potential £j + (iwn) 



in Eq. (1.23) by the effective medium self-energy S^a-/ (i'-^n)- Finally, from comparing 
Eqs. (1.26) and (1.27) we arrive at the desired self-consistency condition 



1 



1 



Ns ^ IUJn - £fc + M - ^AM («Wn) l^n + fJ- - A (^a;„) - T,AAI (jw„) 



(1.28) 



which would give an improved hybridization function A (ia;„) in an iterative procedure. 
We should note in passing that the average medium self-energy and the average Green's 



function (1.25) are the key ingredients in the calculation of the conductivity, which 



as mentioned before involves no vertex corrections within DMFT (Dobrosavljevic and 



Kotliar, 1994) 



It is important to know what the limitations of this approach are. The main one 
is its inability to describe the disorder-induced Anderson metal-insulator transition 



(Anderson, 19581. As the self-consistency condition makes quite apparent, the central 



order parameter of this mean field theory is the average local Green's function, see 
Eqs. (1.21), (1.261 or (1.28). However, as explained by Anderson in the original 1958 



paper (Anderson, 19581, the average local Green's function, which is non-critical and 



finite at the mobility edge, is unable to signal the phase transition between extended 
and localized states. Indeed, the spatial fluctuations of the local Green's function are 
so large that its typical value is far removed from the average one (see Fig. |1.1[ 3 and c) . 
Thus, DMFT cannot describe the Anderson localization transition and one needs to 
go beyond this approximation if Anderson localization effects are to be incorporated. 
It had long been known that CPA has no Anderson transition, so this should not come 
as a big surprise. We will show below, however, that one can in fact address the effects 
of localization, while at the same time retaining the local description of all correlation 
effects. 

1.1.5 Applications of the disordered DMFT 

Early applications of the disordered DMFT scheme focused on the phase diagram of 
the disordered Hubbard model. In particular, the fate of the antiferromagnetic phase 



of the clean model when disorder is introduced was investigated (Ulmke et al., 1995 



Singh et al., 1998). Appropriate incorporation of broken symmetry phases, like antifer- 



romagnetism, requires generalizing the procedure of Section |1.1.4| through the intro- 
duction of sub-lattice structure and spin-dependent single particle quantities, which is 



quite straightforward and will not be discussed here (Georges et al., 1996). The main 



findings are a surprising enhancement of the ordering tendencies at weak disorder and 
strong interactions, which was attributed to a peculiar disorder-induced delocalization 
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effect (Ulmke et at, 1995 Singh et al., 1998). More recently, the phase diagram of 



the paramagnetic Hubbard model (suitable for systems with a high degree of mag- 
netic frustration) has been determined, showing a gradual suppression of the region 



of coexistence of metallic and insulating phases found in the clean case (Aguiar et al., 



2005) 



Kondo disorder. A very attractive feature of the disordered DMFT approach is its 
ability to provide full distributions of local quantities, which in turn may have profound 
effects on the low temperature behavior of physical systems. Consider for example the 



ensemble of effective actions in Eq. ( 1.19 ), which, as explained before, can be viewed as 
an ensemble of Anderson single-impurity problems, with the same conduction electron 
bath (Eq. ( 1.6 )), but different impurity-site energies Sj . As is well known, at sufficiently 
strong coupling U, a local magnetic moment can be stabilized at these impurity sites 
at high temperatures (Anderson, 1961). However, below an energy scale set by the 
Kondo temperature T^j, the moments are "quenched" by the conduction electrons 



and form a singlet bound state (or Kondo resonance) ( 


Anderson, 1961 Kondo, 1964 


Yuval and Anderson, 1970 Anderson and Yuval, 1969 


Nozieres, 1974 Wilson, 1975 



Hewson, 1993). The dependence of Tkj on Sj is given by 

Tkj « £'exp(- 



'l/ppJj 



(1.29) 



where D and pp are the conduction electron half band- width and density of states at 
the Fermi level, respectively, and Jj is the local Kondo exchange coupling constant 



(Schrieffer and Wolff, 1966) 



1 



1 



(1.30) 



Therefore, because of the strong exponential dependence of the Kondo temperature on 
the local parameters, a distribution of site energies can give rise to a wide distribution 



of Kondo temperatures (Dobrosavljevic et al, 1992). As a consequence, depending on 



whether a specific site has T^j < T or T^j > T, it will behave as a free spin in the 
former case or as a quenched inert impurity in the latter one, with significant effects 



on thermodynamic and transport properties, see Fig. (1.3). 

This scenario has received strong experimental support in the context of disordered 
heavy fermion systems. This was initially sparked by NMR experiments done on the 
Kondo alloy UCu4_xPdx, whose broad temperature-dependent line- widths were an- 
alyzed in terms of a distribution of Kondo temperatures (Bernal et al., 1995). The 



same distribution was then used to calculate the magnetic susceptibility and the spe- 



cific heat, with very good agreement with the observed behavior (Bernal et al., 1995). 
In this context, this phenomenology has been dubbed the Kondo disorder model. This 
was particularly striking because this system was among many intensively studied 



heavy fermion compounds (Stewart, 2001 Stewart, 20061 whose properties are in ap- 



parent contradiction with Landau's theory of Fermi liquids (Landau, 19576 Landau, 



1957a Landau, 1959). In particular, the magnetic susceptibility showed an approxi- 
mately logarithmic divergence with lowering temperatures, in contrast with the usual 
saturation to a constant found in weakly, or even some strongly correlated Fermi liq- 
uid metals. The reason for the observed anomalous behavior was quite clear within 
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Fig. 1.3 Distribution of Kondo temperatures for the alloys UCu4Pd (long dashed line) 
and UCua.sPdi.s (short dashed line). Spins with Tk in the hatched area (T > Tk) behave 



as effectively free and lead to a singular thermodynamic response. From reference (Miranda 



et al, 19961 



the Kondo disorder model. Indeed, the distribution of Kondo temperatures needed to 
explain the NMR hne-widths was so broad that P {Tk) ~ Pq = const, when Tk < A, 
where A is some low energy scale of the distribution. In this case, no matter how low 
the temperature is, there are always a few unquenched spins left over with Tk < T 



whose contribution to the susceptibility is Curie-like and large (Fig. (1.3)) 



X{T) 



1 

T' 



(1.31) 



1975) 



Thus, by using a fairly accurate parametrization of the Kondo susceptibility (Wilson, 

1 



\Kondo 

(T) 



(1.32) 



T + aTK 

one can immediately find that the bulk susceptibility obtained from an average over 
the local contributions calculated with the empirical P (Tk) distribution is dominated 
by the Iow-Tr- spins (with Tk T ^ A) and is logarithmically divergent 



xiT) 



P{Tk) 
T + aTK 



dT 



K 



^" dTK^lni^ 
T + aTK \T 



(1.33) 



where Tq is a distribution-dependent constant. 

Clearly, the Kondo disorder model found a natural setting within the DMFT ap- 
proach to disordered systems, which put the phenomenology obtained from NMR on 
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a firmer basis (Miranda et ai, 1996 Miranda et al., 1997a Miranda et al., 1997&I. 
Because heavy fermion systems are characterized by a lattice of ions with incomplete 
f-shells, the most appropriate model Hamiltonian is a disordered Anderson lattice 
Hamiltonian 



HAnd 



J2{Ef,-f,)flf,, 



(1.34) 



in usual notation, and in which we have in general assumed that both f- and c-site 
energies, Efj and e^, as well as the local hybridizations Vj between them are random 
quantities, each with its own independent distribution. The effective local action in 
this case reads 



Seff U) = / E ["l (^) (9r + e, ~ m) c,, (r) + (r) (9. + Ef, - pi) (r) 

+ dr f dr' J2 Ul (r) A (r - r') c,. (r')" 
Jo Jo ^ ^ 

Sl (^) (^) + Ufj^ (r) /,t (r) fl (r) (r) 



and the self-consistency condition is analogous to the one in Eq. (1.24) 
where the local c-electron Green's function GJ^ (iw„) is 



fj 



- A {iujn) , 



(1.35) 



(1.36) 



(1.37) 



and the averaging procedure is performed over the random quantities Sj, Vj and Efj. 
Thus, Tft- fluctuations can have several origins in general, as the local Kondo temper- 
ature is affected by £j, Efj and Vj. Within DMFT, it was possible to better justify 
the ad hoc assumptions of the Kondo disorder model. In particular, one could quantify 
the validity of and thus justify the approximation of calculating the bulk susceptibility 
as an average over single-site contributions ( ^Miranda et al., 1996 ). In addition, good 
agreement was also found with the dynamic magnetic susceptibility obtained through 
neutron scattering experiments (Aronson et al., 1995'). Furthermore, going well be- 
yond the simple Kondo disorder phenomenology, the DMFT approach is able to give 



direct information about transport properties. It was found that (Miranda et al., 1996 



Miranda et al, 1997a Miranda et al, 1997&) 



There is a strong interaction-induced renormalization of the disorder seen by the 
conduction electrons. 
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This, in turn, leads to a rapid suppression of the low-temperature Fermi liquid co- 
herence characteristic of clean heavy fermion materials, as a function of increasing 
disorder. 

Finally, when the quasi-particle coherence is completely destroyed and the distri- 
bution of Kondo temperatures develops a finite intercept in the limit of Tr- — )■ 0, 
the non-Fermi liquid thermodynamics described above is accompanied by a non- 
Fermi liquid linear in T resistivity 



p{T)^Po-AT, 



(1.38) 



where A > 0. Like in the case of the thermodynamic properties, the anomalous 
resistivity is also due to left-over low-T^ free spins off which the conduction 
electrons scatter incoherently. 

It was possible to verify that the self-consistency does not lead to a large disorder 
dependence of the hybridization function A (r). As a result, the distribution of Kondo 
temperatures is fairly sensitive to the bare distribution of random parameters Sj, Efj 
and Vj. Going beyond DMFT, as we will discuss later, one finds that this is an artifact 
of the approximations and, in general, self-consistency leads to a much more robust 
dependence on disorder. 

Elastic and inelastic scattering in the disordered Hubbard model. The interplay be- 
tween local correlation effects and transport is a striking feature which is made almost 
obvious by the DMFT scheme. This has been demonstrated in studies of the disor- 



dered Hubbard model, Eq. (1.18), in references (Tanaskovic et al, 2003) and (Aguiar 



et al, 2004), as we now describe. 



The first study (Tanaskovic et al., 2003) was confined to T = and thus addresses 



only the effects of elastic scattering. This was done using the Kotliar-Ruckenstein slave 



boson mean field theory as the impurity solver (Kotliar and Ruckenstein, 1986) (see 
Section 1.2.2 for further details). The relevant question is how interactions renormalize 
the scattering of quasiparticles by the disorder potential at the Fermi level (hence at 
T — 0). In Hartree-Fock theory, the renormalized disorder potential is determined by 
the self-consistently determined distribution of the electronic charge. This, in turn, is 
governed by the charge compressibility if the charge can be assumed to readjust itself 
to the disorder potential in a fashion dictated by linear response theory. For small U, 
the response is that of a good metal and allows for a flexible adjustment of the charge 
to the bare random potential, leading to a weakened renormalized disorder ("disorder 
screening") (Herbut, 2001 1. For strong interactions close to Mott localization, however, 
the charge compressibility is significantly reduced (the metal becomes increasingly less 
compressible) and Hartree-Fock theory predicts poor disorder screening. 

It was shown in (Tanaskovic et al., 20031 that indeed the efficient disorder screening 



predicted by Hartree-Fock theory is recovered by DMFT at weak interactions. How- 
ever, it was found that another phenomenon intervenes and strong disorder screening 
does occur even as the system approaches the Mott transition. The reason why this 
happens is once again related to the peculiarities of the Kondo effect discussed above. 



Indeed, as several DMFT studies have shown (Georges et al., 1996), the Mott tran- 



sition is signaled by the disappearance of the metallic quasiparticles (Brinkman and 
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Rice, 1970). The coherent nature of these quasiparticles exists only within a narrow 



energy range around the Fermi level whose width is set by the Kondo temperature of 
the associated single-impurity problem: as U ^ Uc, Tk 0. In the disordered case, 
there is a distribution of Tr-'s, but all of them vanish at the transition. Now, the value 
assumed by the renormalized disorder potential on a given site is set by the position of 
the local Kondo resonances (within DMFT), which are known to be strongly pinned to 



the Fermi level (Hewson, 1993). Therefore, Kondo resonance pinning strongly reduces 



the bare disorder fluctuations and screens the disorder rather effectively, leading to a 
correlation-induced suppression of the renormalized disorder, in sharp contrast to the 
simple Hartree-Fock prediction. 



Kondo physics again comes in when one looks at inelastic scattering ( Aguiar et al., 
2004 ) . This was done by means of iterative perturbation theory ( Georges and Kotliar, 



1992j ^Zhang et al., 1993) [ Kajueter and Kotliar, 1996| ) (see Section |1.3.2| for more 



details). Indeed, Tk governs also the temperature above which inelastic scattering 
dominates over elastic scattering. It is found that a gradual and mild temperature de- 
pendence of the resistivity is observed in the weakly correlated regime U W <^ D, 
which is reasonably captured by Hartree-Fock theory. As interactions become of the 
order of (or larger than) the Fermi energy U W ^ D, however, a sharper tempera- 
ture dependence sets in. This is due to the strong suppression of the low-temperature 
scales of the associated Kondo impurity problems, which is not well described within 
Hartree-Fock theory. 

Furthermore, varying the disorder strength W at strong interactions {U ^ D) 
leads to vastly different temperature dependences of the resistivity. Here, the distribu- 
tion of Kondo temperatures defines the range over which inelastic processes become 
progressively more dominant. When U ~ W, a wide distribution of Kondo temper- 
atures is generated and a rather slow growth of the resistivity with temperature is 
found. In contrast, in the cleaner case W U, there is a much narrower distribution 
of Tk 's resulting in a sharp temperature dependence of the resistivity, typical of the 
onset of coherence in heavy fermion materials ( Stewart, 1984[ ) (see Fig. (1.4)). 



1.2 Mott-Anderson transitions: Typical Medium Theory 

Although the dynamical mean field theory described above is able to capture many 
features which are expected to be quite independent of its underlying assumptions 
(e. g., a distribution of local energy scales governing both thermodynamic and trans- 
port properties), it is clear its limitations call for improvements at several points. In 
particular, it would be highly desirable to incorporate Anderson localization effects. 
As we have seen, these are conspicuously absent in the original DMFT, as the latter 
is essentially a mean field theory whose order parameter is A (w) , which in turn is 
determined by the average local Green's function, see Eq. (1.24) or ( 1.27||1.28 ). As 
Anderson localization originates precisely in the spatial fluctuations of this quantity, 
this is not enough. Physically, A (w) represents the available electronic states to which 
an electron can "jump" on its way out of a given lattice site. From Fermi's golden rule, 
the transition ("escape") rate to a neighboring site is proportional to the imaginary 
part of A(w). If this is zero at the Fermi energy, the electrons cannot hop out and 
are effectively localized. In a clean system, in which this quantity is the same at every 
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Fig. 1.4 Scattering rate 1/r as a function of temperature for tlie disordered Hubbard 
model witliin DMFT at (/ = 2D = 2Ep for weak {W = 0.5-D, top panel) and strong disorder 
{W — 2D, bottom panel). Insets show the dependence of the local Kondo temperature ~ Zi 
on the site energy Si. At weak disorder, a sharp distribution of Kondo temperatures gives rise 
to a steep rise, whereas stronger disorder leads to a broad distribution and a slow temperature 
dependence. From reference (Aguiar et al, 20041. 

site, this is a good order parameter for localization, as in the case of the clean Mott 
transition. In a highly disordered system, A (uj) shows strong spatial fluctuations from 
site to site and its average value is not a good measure of the conducting properties. A 
"typical" site in an Anderson insulator will have a hybridization function (cj) with 
large gaps and a few isolated peaks, reflecting the nearby localized wave functions 
that have an overlap with it. The vanishing of its imaginary part signals the electron's 
inability to leave the site and is a good indicator of localized behavior. However, aver- 
aging over the whole sample washes out these gaps hiding the true insulating behavior. 
The discrepancy between the typical and average values of A^ (w) persists even on the 
metallic side, where Atyp (uj) can be much smaller than A (w). 

Two alternative routes can be taken at this point. The ideal solution is to track the 
actual local hybridization or escape rate at each site. We will focus on this possibility in 
Section |1.3[ where we analyze the so-called Statistical Dynamical Mean Field Theory. 
The other option is to focus on a simpler, yet meaningful measure of the escape rate 
which, though incapable of incorporating the richness of the actual local realizations. 
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does not "throw away the (locahzation) baby with the bath water". Here we should 



take as guidance the remark by Anderson that "no real atom is an average atom" ( An- 



derson, 1978 1 and seek a more apt description of a "real atom' 
of the typical escape rate ImAt^p (w 



Indeed a good measure 
is the geometric average exp (In [ImA^ (w)]). This 



has the great advantage of serving as an order parameter for the localization transi- 
tion: indeed, in contrast to the algebraic average, the geometric average vanishes at 



the mobility edge (Anderson, 19581 (see also (Janssen, 1998 Mirlin, 2000 Schubert 



et ai, 2010)). We can then reason by analogy with the regular Dynamical Mean Field 



Theory approach explained above and construct a self-consistent extension centered 
around this typical escape rate function. This theory has been dubbed the Typical 
Medium Theory (TMT) ( [Dobrosavljevic et ai, 2003a[ ). 



1.2.1 Formulation of the theory 



We can proceed by analogy with the DMFT equations as explained in Section 1.1.4 
see Eqs. ( 1.25|1.28 l, to obtain the Typical Medium Theory. We again focus on a 



disordered Hubbard model, Eq. (1.18), and imagine replacing the disordered medium 
by an effective typical medium described by a self-energy function 'Stmt i'-^)- How 
do we determine this self-energy? Focusing on a generic site j, it is described by 



an effective action which has the same form as Eq. (1.19). The hybridization function 



A (w) is still left unspecified at this point, but we envisage that it will reflect a "typical" 
site as opposed to an "average" one, so we set A (cj) ~ \yp (i^) in Eq. (1.19). The 



local Green's function is still defined as in Eqs. (1.22) and (1.23). The loca' 



density of 



states is given by the imaginary part of the local Green's function 

pj (oj) = —ImGjj (oj — iS) . (1.39) 

TT 

The typical local density of states can be defined through its geometric average 



ptyp (w) = exp 



dSjP (e j) In Pj {uj) 



(1.40) 



Note that we have reverted to the real frequency axis because we need a positive- 
definite quantity in order to be able to define a geometric average. To preserve causality, 
the typical local Green's function is obtained through the usual Hilbert transform 



Gtyp (w) 



(1.41) 



Note the analogous average quantity in the second equality of Eq. ( 1.26 ), which appears 



in DMFT. The typical medium self-energy Htmt (w) is then defined through the 
inversion of the following equation 



Gtyp (w) 



1 



^TMT 



(1.42) 



which is the analogue of the first equality in Eq. (1.26). Finally, the loop is closed by 
setting 
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(1.43) 



+ /i - Atyp (cj) - Etj\/t (i^) 
which can be used in an iterative scheme to generate an updated hybridization function 



Atyp (w) and is the analogue of Eq. ( 1.28 ). It becomes clear that the crucial difference 



between TMT and DMFT is the replacement of the average local Green's function by 
the typical one (see Fig. |l.l[ 3 and c). This has been shown to capture even quantitative 
features of the Anderson localization transition, as we will discuss below. For a more 
comprehensive review, see ( [Dobrosavljevic, 2010 ). 



1.2.2 Applications of TMT 

We will now describe the most important results obtained from the TMT theory of 
disordered systems. We will focus on the non- interacting case in Section 1.2.2| and on 
the disordered Hubbard model in Section [l. 2. 2[ We should also mention a study of the 



Falicov-Kimball model within TMT in reference (Byczuk, 2005). 



Critical behavior in the non-interacting case. As a first test of the usefulness of the 
TMT approach, it was first applied to the non-interacting three-dimensional case, 
Eq. (1.18) with [7 = (Dobrosavljevic et al., 2003a). In fact, the results of ap- 
plying the TMT equations to a cubic lattice were directly compared to a numer- 
ical diagonalization of the Hamiltonian. In particular, the numerically determined 
arithmetic and geometric averages of the local density of states at the Fermi level, 
Pai) (w = 0) and pgeo (w = 0), were compared to the results of CPA (Elliott et al, 



1974[ Economou, 2006) and of TMT (see Fig. 1.5). A remarkably accurate agreement 
between pav (w = 0) and CPA was observed. As is known, this quantity is not critical 
at the Anderson transition. On the other hand, the geometric average of the local 
density of states does vanish at a critical disorder strength Wc- A reasonably good 
agreement between the numerical pgeo (w — 0) and TMT is found for most values of 
the disorder strength, even though TMT misses the correct critical behavior. This is 
not too surprising as TMT has the flavor of a mean field theory. Physically, it is clear 
the Pgeo (w) should be viewed as the density of extended states of system, decreas- 
ing with increasing disorder and eventually vanishing altogether for sufficiently large 
randomness. Thus, the spectral weight described by pgeo (w) is not conserved. 

In fact, further insight into the critical behavior of TMT can be achieved ana- 
lytically ( Dobrosavljevic et al., 2003a[). By assuming an elliptic density of states for 

the clean lattice, Poiy) = ~ (ij)^' '-^'^ proved that in the critical re- 

gion W — )■ Wc = eD, the typical density of states (given, as usual, by the geometric 
average) assumes a universal form 



where 



Ptyp (w, W)^p{oj^O,W)f [UJ/UJO (W)] , 



p(c^ = 0,W^)= ( - ) {Wc-W), 



(1.44) 
(1.45) 



the frequency scale 
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Fig. 1.5 Comparison between tiie TMT and CPA witii numerical results (circles) on the 
non-interacting 3D tight-binding model with diagonal disorder. The average local density 
of states at the band center pav is remarkably well captured by the Coherent Potential 
Approximation (dashed line) but is not critical at the Anderson localization transition at 
Wc ~ 1.38 (in units of the clean bandwidth). The geometrically averaged local density of 
states ptyp, on the other hand, vanishes at the transition and is well described by the Typical 
Medium Theory (full line), except for the detailed critical behavior. The inset shows the 
numerics for the critical behavior of ptyp as a function of the correlation length ^ (W), where 
C (W) ~ I W - Wol", V ^ 1.58. From reference pobrosavljevic et al, 2003a[). 



W), 



(1.46) 



1 



Note that Eq. (1.451 gives 



and the scaling function has a simple form / (x) 
an order parameter critical exponent (3tmt = 1, which should be compared to the 
accepted value in three dimensions « 1.58 (Slevin and Ohtsuki, 1999). 



The disordered Hubbard model at half hlling. We now direct our attention to disor- 
dered interacting systems. The TMT was applied to the disordered Hubbard model 
(Byczuk et al., 2005 Aguiar et al, 2009). The phase diagram of the paramagnetic 



half-filled case at T = was obtained with two different impurity solvers: Wilson's 



numerical renormalization group (NRG) (Wilson, 1975) was used in (Byczuk et al, 



2005) and the Kotliar-Ruckenstein slave boson mean field theory (SB4) (Kotliar and 



Ruckenstein, 1986) was applied in (Aguiar et al., 2009). The results obtained largely 
agree with each order but there are some small discrepancies in the details. Essentially, 
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three different phases are observed: a disordered correlated metal phase, characterized 
by pc = Ptyp (t^ = 0) 7^ 0, a Mott-like insulating phase (which we will call simply a 
Mott insulator) and a Anderson-like insulating phase (for which we will use the name 
Mott- Anderson insulator), the latter two phases having = 0. We will discuss each 
set of results and then the discrepancies. For reviews, see (Dobrosavljevic, 2010) and 
( Byczuk et ai, 2009 6| ). 

For small values of disorder and interaction, both approaches find that the system 
is metallic, although this is not too surprising. There is also agreement on the fact that, 
for a fixed small disorder strength, the order parameter pc increases with increasing 
interaction {U > W). This is a reflection of disorder screening by interactions (see 
Section 1.1.5). Eventually, at a critical value of the interaction strength Uc{W), the 
order parameter pc exhibits a finite jump and drops to zero, signaling a metal-Mott 
insulator phase transition. Furthermore, both methods agree that, for a fixed small 
value of U, as one increases the disorder {W > U) pc decreases monotonically, indicat- 
ing that the spectral weight due to extended states is decreasing. At the critical value 
Wc {U), Pc vanishes and the system enters a Mott- Anderson insulating phase. 

The differences in the results of the two approaches are the following. The NRG- 
based TMT (Byczuk et ai, 2005) predicts the metal-Mott insulator transition for U > 
W to be first order in character, with typical hysteretic behavior: the metallic solution 
is locally stable for U < Uc2 {W) and the Mott-insulating solution is locally stable for 
U > Uci [W), where Uci {W) < Uc2 {W). In the coexistence region Uci [W) < U < 
Uc2 {W), both solutions can be stabilized (although only one is a true global energy 
minimum at each U). The SB4-based TMT ( [Aguiar et al, 2009| ), on the other hand, 
predicts no hysteresis. To understand this, it should be mentioned that SB4 is a good 
description of the low-energy properties of the impurity spectral function, although it 
misses the higher-energy features which give rise to the Hubbard bands in the lattice. 
An important ingredient of this description is the quasiparticle weight = Z (si), 
well known from Fermi liquid theory, which in the impurity problem - Eq. (1.19) 
- determines the width of the Kondo resonance (essentially the Kondo temperature) . 
Formally, it appears in the local Green's function at site i 

lUJn - e,+ p - Z^Atyp (lUJn) 

where ii is the renormalized local site energy, which gives the position of the resonance. 
In the clean lattice case, Z vanishes continuously as the interaction strength is tuned 
to its critical value U — >■ Uc, Z ^ Uc — U — >■ 0, signaling the Mott localization of 
the itinerant electrons, which become localized magnetic moments. The continuous 
nature of the transition in the clean case, with no accompanying hysteresis, survives 
the introduction of disorder within the SB4-based TMT. In fact, in this approach all 
the Zi^s vanish at a unique Uc{W). It should be stressed that in both approaches 
Pc is discontinuous at the transition, a fact which can be ascribed (at least for small 
disorder) to the observed perfect screening of disorder {pav (w = 0) — >■ pgco (w = 0) as 
U ^ Uc {W)) and the pinning of the clean density of states at the Fermi level to its 
non-interacting value (see Section 1.1.5). 

Moreover, the results obtained with the NRG impurity solver indicate that the 
transition in the region W > U is such that Pc ^ continuously as W ^ Wc (U) 
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(Byczuk et al., 2005). This is in contrast to the SB4-based approach, which finds that 
Pc exhibits a discontinuous jump to zero at Wc{U) (Aguiar et al., 20091. In fact, 
the latter method brings out an important ingredient which significantly enhances 
the physical understanding of the TMT approach to this problem. This is achieved 
by tracking the behavior of the quasiparticle weights Zi. For a given fully-converged 
hybridization function Atyp («Wn), the ensemble of impurity problems is characterized 
by the function Z {si) for \ei\ < W/2 (a uniform disorder distribution is assumed). 
This function has the property that, as the phase transition is approached, Z (ei) —, 
for \ei\ < U/2, whereas Z (e^) -)> 1 if |g^| > U/2 (see Fig. [Tel ( [Aguiar et~. 
2006). In the SB4 language (Kotliar and Ruckenstein, 1986), Z — >■ implies a singly 
occupied site with a localized magnetic moment, whereas Z — >■ 1 means either a doubly 
occupied or a singly occupied site, either of which is essentially non-interacting. Thus, 
in the region W > U, a fraction of the sites, those with \ei\ < U/2, experience Mott 
localization, while those sites with \ei\ > U/2 undergo Anderson localization. The 
picture that emerges is that of a spatially inhomogeneous system, composed of Mott- 
localized droplets intermingled with Anderson insulating regions. This situation has 
been dubbed a "site-selective Mott transition" (Aguiar et al., 2009). Analytical insight 
into the SB4-results can be brought to bear in order to show that in that approach 
any finite U renders the vanishing of pc discontinuous, in sharp contrast to the non- 



interacting case (Dobrosavljevic et al., 2003a 



Finally, in the intermediate region oiW k,U , it was suggested, based on the NRG 
results that there might be a crossover from a metal to a disordered Mott insulator 
(Byczuk et al., 2005). However, since their results clearly show regions where Pc 7^ 
and regions where Pc = 0, we believe the correct interpretation is to identify the 
former as metallic and the latter as insulating. Besides, this is expected from the 
sharp distinction between an insulator and a metal at zero temperature. Having said 
that, however, it is clear that the nature of the transition in this problematic region 
certainly deserves a further more detailed investigation. 

The disordered Hubbard model was also studied within TMT by allowing anti- 
ferromagnetic order on a bipartite lattice (Byczuk et al., 2009a). The phase diagram 
at zero temperature as a function of disorder and interactions was determined, with 
the identification of both paramagnetic and antiferromagnetic metallic phases (for fi- 
nite disorder only), an antiferromagnetic Mott insulating phase and a paramagnetic 
Anderson insulating phase at large disorder strength. 



1.3 Mott-Anderson transitions: Statistical DMFT 

As outlined above, the most natural and accurate extension of the DMFT philosophy 
which can incorporate Anderson localization effects is one which replaces the algebraic 
average of the hybridization function A (w) (as in the original DMFT) or its typical 
value/geometric average ^typ (w) (as in the TMT), by the actual realizations of Aj (w) 
at each site ( Dobrosavljevic and Kotliar, 1997[ Dobrosavljevic and Kotliar, 1998) (see 
Fig. |1.1[ 3 and c). As is to be expected, the complexity of the equations increases 
considerably and one has to rely heavily on numerical computations. However, many 
insights have been obtained regarding the systems analyzed. Besides, a much larger 
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Fig. 1.6 Z (si) showing the two-fluid behavior of the Mott-Anderson transition within 
TMT. As the transition is approached (t 0) for W = 2.8 > i7 = 1.75, sites with 
\ei\ < (7/2 ~ 0.31W represent Mott insulating regions with Z (si) — > 0, whereas sites with 



\ei\ > U/2 are Anderson localized with Z (si) -> 1. From reference ( Aguiar et ai, 2006| 



degree of universality of the distributions is observed when compared with the much 
more "rigid" approaches of DMFT or TMT. 

1.3.1 Formulation of the theory 

Let us examine how the statDMFT works. We begin with the by now famihar disor- 



dered Hubbard model of Eq. (1.18) and focus, as usual, on the dynamics of a given 



site J, dictated by an effective action (we now revert back to imaginary time) 

„ Jo Jo 



rP 

+ U I drnj^ (t) n-n (r) 
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Notice the crucial difference now: the hybridization function Aj (iiVn) is now site- 
dependent. Each site, besides having a different energy Sj, also "sees" a different local 
environment Aj (iujn)- The local dynamics is again encoded in a site-dependent self- 



energy T,j (iwn), obtained from the local Green's function as before, see Eqs. (1.221 



and 1.23). It is important to note that this description assumes a diagonal (albeit site- 



dependent) self-energy function, adhering to the generic philosophy of incorporating 
only local interaction effects. 

Unlike the previous DMFT or TMT approaches, statDMFT does not try to mimic 
this self-energy function {iuJn) through any type of effective medium. Instead, we 
choose to include its full spatial fluctuations. We do so by appealing to the physical 
picture of the self-energy as a shift of the local site energy, albeit a complex, frequency- 
dependent one: ej — >■ Sj + (?cij„). Single-particle propagation can thus be viewed as 
described by the effective resolvent 



G (iidn) = il-^n^ ^ t — £ — E (iuJn) 



(1.49) 



where 1 is the identity operator, t and e are, respectively, the hopping and site-energy 



terms (first and second terms of the Hamiltonian (1.18)), and the matrix elements of 
the self-energy operator E (iiOn) are given in the site basis as 



E (ZW„) i) {iuJn)Sr 



(1.50) 



Although any matrix element of the resolvent (1.49), both intra- and inter-site, can in 



principle be calculated (which is important, for example, for a Landauer-type calcula- 
tion of the conductivity), the self-consistency requires only the diagonal part, related 
to the local Green's function 



G{iu}n) j) = Gjj (iujn) = 



1 



- £j - ^3 



Ej (iuJn) ' 



(1.51) 



This last equation closes the self-consistency loop by providing, in an iterative scheme, 
an updated hybridization function for each site Aj {iLu„). We summarize the self- 
consistency loop for completeness: 

1. For a given realization of disorder, Sj, start from a set of "initial trial" hybridiza- 
tion functions A^- (ia;„), one for each site. 

2. For the set of effective actions ( 1.48 ), calculate the local Green's function Gjj (iujn) 
(Eq. (1.22)) and the local self-energy Ej (za;„) (Eq. (1.23)) for each site. 

3. Invert the matrix resolvent (1.49) and get its diagonal elements Gjj (iaj„). 

by equating these di- 



Obtain an updated set of hybridization functions Aj {iuj„ 
agonal elements to the expression of the local Green's functions, Eq. (1.511 



In general, the set of Eqs. ( 1.48|1.5T | forms the so-called Statistical Dynamical Mean 
Field Theory of disordered correlated electron systems. It should be noted that, besides 



the challenge of solving the ensemble of impurity problems represented by Eq. (1.48) 



Eq. (1.51) poses the numerical problem of inversion of a complex matrix for each 



value of the frequency, which can be very time-consuming. The pay-off is a description 
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which incorporates all Anderson localization effects. Indeed, when interactions are 



turned off, the theory becomes exact, since, e. g., Eq. (1.49) becomes the exact single- 
particle Green's function (from which transport properties can be obtained with the 
Landauer formalism). In the absence of randomness, we recover, of course, the DMFT 
equations. In the presence of both disorder and interactions, this is the optimal theory 
of disordered interacting lattice fermions which includes only local correlation effects. 

1.3.2 Early implementations of statDMFT: the Bethe lattice 

It had long been known that the Bethe lattice (or "Cayley tree") leads to consider- 
able simplifications of the treatment of non-interacting disordered systems. For ex- 
ample, the so-called self-consistent theory of localization of Abou-Chacra, Anderson 



and Thouless (Abou-Chacra et al, 1973) becomes exact on a Bethe lattice. This is 



so because the local Green's function with a neighboring site removed, Eq. (1.4 1, sat- 
isfies a single compact stochastic equation in that lattice, which allows for quite an 
efficient analysis, both analytically and numerically. It was shown, for example, that 
for a coordination number z > 2, there is always an Anderson transition at a non- 
zero critical disorder strength Wc- This transition has been extensively studied fMirlin] 



and Fyodorov, 1991 ) and is regarded as a large-dimensionality limit of the Anderson 
transition. Although it has an anomalous critical behavior, with an exponential rather 
than a power-law dependence, this description has been very fruitful, specially if one 
is interested in the non-critical region. 

Given the complexity of the full statDMFT equations, this suggested that the pre- 
liminary investigations could be carried out on a Bethe lattice. Here, it is appropriate 
to comment on what is, in our view, a misunderstanding of the conceptual basis of the 



statDMFT approach. It has been stated (Semmler et al., 20106) that there is some- 



how a conceptual difference between the statDMFT as applied to the Bethe lattice 
and the statDMFT used to analyze realistic lattices, like the square or cubic ones 



(for these, see the later Section 1.3.3). The misunderstanding comes from assuming 
that whereas on a realistic lattice one has a fixed disorder realization, which is then 
solved by statDMFT, therefore defining a deterministic problem, on the infinite Bethe 
lattice one does not deal with fixed disorder realizations but with distributions, and 
the approach becomes non- deterministic and "statistical" (we realize the origin of the 
misunderstanding may have been the use of this word in the name of the method 



(Dobrosavljevic and Kotliar, 19981). This distinction is unfounded if one realizes that 
on any infinite lattice with a single fixed disorder realization (with random, spatially 
uncorrelated, parameters) , the infinite values of local quantities (such as the local den- 
sity of states) on each lattice site give rise to a statistical distribution of local quantities 



(assuming the lattice-translational invariance of the distributions, see (Lifshits et al., 
[1988) for a careful discussion). Thus, when one solves a given Hamiltonian with stat- 
DMFT on an infinite Bethe lattice, one is actually solving, in practice as well as in 
principle, for a single fixed disorder realization. Each iteration of the method outlined 



in (Abou-Chacra et al, 1973) corresponds to "going outwards" on the branches of 
a fixed disorder realization of an infinite Bethe lattice, while at the same time accu- 
mulating random variables and building a histogram of local quantities. By the same 
token, if one could solve the statDMFT equations for a single disorder realization on 
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an infinite realistic (say, square) lattice, each one of its sites would contribute one 
random variable for a histogram of the same local quantities (which would be, obvi- 
ously, different from the ones obtained from the differently connected Bethe lattice). 
In practice, of course, one solves many disorder realizations of finite, hopefully large, 
realistic lattices in order to generate distributions with good statistics. However, fun- 
damentally, there is no conceptual difference between the statDMFT solutions on the 
two types of lattices. 



The Mott-Anderson transition. In (Dobrosavljevic and Kotliar, 1997 Dobrosavljevic] 



[and Kotliar, 1998) , the first implementation of the statDMFT theory, as applied to 
the Mott-Anderson transition described by the disordered Hubbard model was com- 
pleted. This was done by using a Fermi liquid parametrization of the associated zero- 
temperature impurity problems, namely, the infinite-C/ slave boson mean-field theory 



(Read and Newns, 1983 Coleman, 19871. Like its Kotliar-Ruckenstein finite-C/ coun- 
terpart (Kotliar and Ruckenstein, 1986), this theory captures the low energy sector 
and is known to give a quantitatively good description of this limit. This Fermi liquid 
description is encapsulated in just two parameters: the quasiparticle weight Zj and 
the effective level energy (or Kondo resonance location) Sj. The local self-energy is 
written as 

(iujn) = (1 - Zr^) iLo„ + ^ - £j + M, (1.52) 
leading to a local Green's function 



The results of (Dobrosavljevic and Kotliar, 1997 Dobrosavljevic and Kotliar, 1998) 



reveal that, as in the non-interacting case, the Anderson-Mott transition can be iden- 
tified by the vanishing of the typical local density of states (as described by the ge- 
ometric average) at a certain critical disorder strength. Interestingly, in contrast to 
non-interacting electrons, the critical behavior is conventional (power-law) and the 
average density of states is divergent. There is at present no good understanding of 
this divergence. 

Furthermore, a great opportunity afforded by statDMFT is the ability to inves- 
tigate distributions of local quantities. By tracking the distribution of quasiparticle 
weights, it was found that it broadens considerably with increasing disorder, showing 
a characteristic power-law form at large randomness 

P(Z)-Z"-\ (1.54) 

with the exponent a = a (W) a smooth function of disorder. It should be remembered 
that Zj determines the local Kondo temperature Tj^j ^ Zj and thus governs the 
local contribution to thermodynamic quantities such as the magnetic susceptibility 



and specific heat (see Section 1.1.5). In fact, by averaging over this distribution of Z's 
as in Eq. ( 1.33[ ), one finds power-law dependences for these quantities as well 



x(r)~C(r)/T~T"~i. (1.55) 
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As disorder increases, a decreases and eventually becomes smaller than 1 well be- 
fore the Mott- Anderson transition. When this happens, the thermodynamic response 
becomes singular and non-Fermi liquid-like 



X(T^ 0) ^ oo. 



(1.56) 



Many different correlated systems have indeed been shown to exhibit this form of 
anomalous behavior (Stewart, 2001), with non-universal exponents a. The presence 



of non-universal, smoothly varying exponents characterizing divergences in physical 
quantities is reminiscent of a large class of disordered systems and is usually dubbed a 



quantum Griffiths phase (for reviews, see ( |Miranda and Dobrosavljevic, 2005 Vojta, 
2006)), by analogy with a similar situation in classical systems first analyzed by Grif- 
fiths (Griffiths, 1969). Most other known examples of quantum Griffiths phases had 



been found in the vicinity of magnetic phase transitions in the presence of disorder. 



most notably in insulating magnets (Fisher, 1992 


Fisher, 1995 


Guo et al, 1996 Pich 


et al, 1998 


Motrunich et al, 2001 


1, but also in metallic systems ( 


Castro Neto et al, 


1998 


de Andrade et al, 1998 


Castro Neto and Jones, 2000,,Millis et al, 2002 1. Here, 



however, the characteristic power laws are found in the vicinity of the paramagnetic 
Anderson metal-insulator transition and hence the name Electronic Griffiths phase 
was adopted. This anomalous behavior is apparently not at all dependent on the par- 
ticular details of the disordered Hubbard model. Very similar power-law distributions 
of Kondo temperatures were also found in Bcthe lattice implementations of statDMFT 



for the disordered Anderson lattice Hamiltonian (1.34 1 (see Fig. ( 


1.7 


I) ( 


Miranda and 


Dobrosavljevic, 1999 


Miranda and Dobrosavljevic, 2001 


Aguiar et al, 2003 





It should be remembered that the forms of the distributions of Kondo temperatures 
obtained within DMFT were strongly dependent on the shape of the bare distributions 
of parameters. This is in sharp contrast to the ubiquitous power-law distributions found 
in the statDMFT approach. Thus, although the exponent a is disorder-dependent and 
non-universal, the power-law shape is quite independent of whether the bare parame- 



ters are given by, say, uniform, Gaussian or binary distributions (Aguiar et al, 2003). 
This is again easy to understand if we note that the local Kondo temperature depends 
exponentially on the local density of states at the Fermi level (and less strongly on 
its value at higher energies). Now, due to the extended nature of the electronic wave 
functions in metallic systems, the density of states at one site is influenced by spatial 
fluctuations at very distant sites and thus samples a great number of local environ- 
ments. The resulting distributions of local quantities thus reflect this long-distance 
sampling. 

In order to understand why this effect leads specifically to a power law, an effective 
model was proposed in (Tanaskovic et al, 2004). The effective model consisted of 



a disordered Anderson lattice model with Gaussian distributed conduction electron 
disorder {sj in Eq. (1.34)) 

1 



exp I 



]/2W^) , 



(1.57) 



'2tiW 

treated within DMFT. In DMFT, the hybridization function A (iw„) is site-indepen- 
dent and the Kondo temperature distribution is solely determined by fluctuations of 
It can then be shown that (Tanaskovic et al, 2004) 



Mott-Anderson transitions: Statistical DI\/IFT 29 



1.0 




I ' 1 ' <— 

-12.5 -7.5 -2.5 

x=log(TJ 



Fig. 1.7 Power-law distributions of Kondo temperatures in a disordered Anderson lattice 
model for different levels of disorder: the linear behavior for small values of log Tk implies a 
power-law P [Tk) ~ T^~^. The inset shows the exponent a as a function of disoder W. From 



reference (Miranda and Dobrosavljevic, 20011 



(1- 



and A is determined by other model 



where is the Kondo temperature at Sj 
parameters but does not depend on Sj. It is easy to show from Eqs. (1.571 and (1.581 
that 



PiTK)=P[e,{TK) 



dTt 



K 



K 



(1.59) 



which is precisely the power-law distribution of Kondo temperatures found generically 
within statDMFT treatments. We note, in passing, that this kind of argument is generic 
to all known quantum Griffiths phases: the relevant energy scales are exponentially 



suppressed by a certain random parameter (see Eq. ( 1.58 )), whose probability is in turn 
also exponentially small (see Eq. (1.57)) (Miranda and Dobrosavljevic, 2005 Vojta, 



2006) 



What is the relation between these results, obtained within DMFT, and the power 
laws observed in the applications of statDMFT? In statDMFT, the single (average) 
hybridization function of DMFT gets replaced by a strongly fluctuating distribution of 
local hybridizations Aj (ia;„). The imaginary part of each of these functions describes 
the available density of states for Kondo screening at site j and enters the expression 
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for the local Kondo temperature much like pp does in Eq. (1.29). From the central- 
limit theorem, fluctuations of the available densities of states around the mean value 
are generically Gaussian for weak and intermediate disorder and lead to a power-law 
distribution of Tr-'s in a fashion quite similar to the effective model. Indeed, detailed 
calculations showed that the effective model is quite accurate when compared with full 
statDMFT results (Tanaskovic et al., 2004). Crucially, these arguments can be used 



to show that the non-Fermi liquid behavior occurs already at quite moderate values of 
disorder and strictly precedes the Anderson metal-insulator transition. This elucidates 
then the microscopic origin of the electronic Griffiths phase. 

Iterative perturbation theory as impurity solver. Most of the early statDMFT results 
on the Bethe lattice were obtained through the use of the infinite-U slave boson mean- 



field theory (Read and Newns, 1983 Goleman, 1987) as impurity solver. These are 



good descriptions of the low-energy coherent Fermi-liquid part of the impurity spec- 
trum but fail to account for inelastic scattering at low energies as well as higher-energy 
incoherent features such as upper and lower Hubbard bands. A technique which is able 



to incorporate there features is the so-called iterative perturbation theory (Georges 



and Kotliar, 1992 Zhang et al., 1993 Kajueter and Kotliar, 1996). The iterative 



perturbation theory approach lends itself also more easily for an analysis of the tem- 
perature dependence of physical quantities. It has been used to analyze the disordered 



Anderson lattice model (Aguiar et al, 2003) as well as the disordered Hubbard model 



(Semmler et al., 2010a). It suffers from the disadvantage of not being able to capture 



the correct exponential dependence of the low-energy Kondo scale, see Eqs. (1.29) and 
( |1.58 ), leaving out, therefore, the possibility of characterizing Griffiths phase behavior 
(Section [TX2| . 

In the case of the disordered Anderson lattice, one of the interesting findings was 
the interplay between elastic scattering off the disorder potential and inelastic electron- 
electron scattering (Aguiar et al., 2003). If one uses the inverse of the typical density of 



states at the Fermi level as a rough guide to the resistivity (a direct calculation of the 
resistivity was numerically prohibitive at the time of those studies), its temperature 
dependence is found to be quite sensitive to the amount of disorder. Indeed, low- 
disorder regimes are marked by an increase of the resistive properties with rising 
temperatures, signaling the onset of inelastic scattering processes, much like in the 
clean case. On the other hand, strongly disordered samples exhibit a decrease in the 
resistivity with increasing temperatures, because a decrease in the effective elastic 
scattering outweighs the increase in the inelastic one. This type of fan-like family of 



resistivity curves (see Fig. 1.8) has been seen in several disordered strongly correlated 
materials, being known as Mooij correlations (Mooij, 1973). They are seen to arise 



here within a local approach to electronic correlations and disorder. 

More recently the Hubbard model with binary alloy disorder has been also studied 



on a Bethe lattice with iterative perturbation theory as impurity solver ( Semmler et al., 



2010a). Particular attention has been paid to the dependence of the distribution of local 



densities of states on the small imaginary part ( "broadening" ) that is usually added to 
the frequency in numerical determinations of Green's functions. The dependence of the 
distribution of local densities of states on this parameter can be used to characterize the 
localized or extended nature of the electronic states. The main result of that paper is 



Mott-Anderson transitions: Statistical DMFT 31 



1.4 



1.2 



^ 1.0 



Q. 



0.8 



0.6 



o- 



-OW= 1.4 



□ Q W = 1 .2 

0---0W= 1.1 
^ — AW = 0.71 
^--<]W = 0.35 
V VW = 0.25 




0.4 



0.000 



0.001 



0.002 



0.003 



T 



Fig. 1.8 Temperature dependence of the inverse typical conduction electron density of states 
of the disordered Anderson lattice model for different levels of disorder. This "fan-like" family 
of curves is generic to many disordered, strongly correlated systems ("Mooij correlations"). 



From reference ( Aguiar et al. , 2003 1 



the determination of the zero-temperature phase diagram of the model at a particular 
filling as a function of interaction and disorder strengths and the identification of 
regions with metallic, Mott-Anderson insulating and band insulating behaviors. In 
particular, the opening of the Mott gap occurs at values of the interaction and disorder 
strengths for which the gapless system is already Anderson localized. It is difficult 
to compare this phase diagram with the one obtained by the infinite-U slave boson 



mean field theory (Dobrosavljevic and Kotliar, 1997 Dobrosavljevic and Kotliar, 1998) 



because the models were solved with different types of disorder and in different regimes. 
1.3.3 StatDMFT on realistic lattices 

Although the Bethe lattice implementations were very informative, its exotic con- 
nectivity introduces some unwanted features, especially with regard to the critical 



behavior of the Anderson transition (Mirlin and Fyodorov, 19911, which is believed 



to correspond to some infinite-dimensional limit. Therefore, it is important to check 
what the effects of finite dimensions are. Motivated by this, the full statDMFT equa- 
tions have been implemented in realistic lattices in recent years. Given the successful 
application of DMFT to the description of the Mott-Hubbard transition, a natural 
candidate is the analysis of the disordered Mott-Hubbard transition. 
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The disordered Mott-Hubbard transition. As with any other phase transition, the 
characterization of the effects of disorder on the Mott transition poses an important 
yet difficuh problem. In the particular case of phase transitions incontrovertibly de- 
scribed by an order parameter, this characterization has seen considerable advances. 
In insulating quantum magnets, many examples have been found in which a whole 
vicinity of the disordered critical point is described as a quantum Griffiths phase. 
This is a phase in which rare regions of nearly ordered material dominate the physics 



and the thermodynamic response becomes divergent (Miranda and Dobrosavljevic, 



2005 Vojta, 2006). The sizes and energy scales governing the rare regions span sev- 



eral orders of magnitude and their description requires taking account of very broad 
distributions. Furthermore, in many cases, as the critical point itself is approached, 
the relative widths of the distributions grow without limit, a situation generically de- 
scribed as an infinite randomness fixed point. This has been well established in systems 
with both Ising and continuous symmetry ( Fisher, 1992| Fisher, 1994 Fisher, 1995 



Hyman et al, 1996; "Yang et al.,l996': 'Guo et al, 1996; Hyman and Yang, 1997[|Fisher1 



and Young, 1998;,Pich et al, 1998, , Narayanan et ai, 19996, .Narayanan et ai, 1999a[ 



Motrunich et al., 2001 Refael et al., 2002). More recently, a symmetry-based clas- 



sification scheme of these Griffiths phases has been proposed and applied to several 



different systems with great success ( 


Vojta, 2003 


Vojta and Schmalian, 20056 


Vojta 


and Schmalian, 2005a Vojta, 2006 |Hoyos and Vojta, 2006 


Hoyos et ai, 2007 


Hoyos 


and Vojta, 2008| Vojta et ai, 2009 


). The Mott transition, however, poses a problem 



of a different nature, as it is not described by an order parameter in an obvious way. 
It is thus not clear how to extend the above insights into its description. 

The implementation of statDMFT offers a natural way out. In particular, the clean 
problem is aptly described by DMFT, as we mentioned at the end of Section 1.1.2 The 
first implementation of statDMFT for the disordered Hubbard model was performed 



in (Song et al., 2008) using the so-called "Hubbard I" (HI) approximation (Hubbard, 



1963 ) as the impurity solver. This simple impurity solver captures the physics close to 



the atomic limit and is therefore convenient for a study of the Mott insulating phase. 
However, it suffers from the deficiency of not giving rise to a quasiparticle peak in 
the single-impurity spectral function, a feature known to exist for any finite U when 
the hybridization function is that of a good metal. The method was applied to the 
cases of half and quarter filling and the density of states and inverse participation ratio 
were obtained. Interestingly, the localization length has a non-monotonic behavior as a 
function of the interaction U: it increases initially with U, showing a tendency to delo- 
calization, but eventually decreases and becomes even smaller than the non-interacting 



value at the Mott transition. At quarter filling, an Altshuler-Aronov (Altshuler and 



Aronov, 1979) density of states anomaly is found, which is however absent at half- 



filling, due to an interaction induced-suppression of the charge susceptibility. 

Shortly afterwards, the disordered Hubbard model in a two-dimensional square 



lattice at T = was solved with statDMFT (Andrade et al., 2009a Andrade et al.. 



20096) using the Kotliar-Ruckenstein slave boson mean field theory as the impurity 



solver (Kotliar and Ruckenstein, 1986). Very similar results were also obtained (Pez- 



zoli et al., 2009 Pezzoli and Becca, 2010) within an approach based on a Gutzwiller 



variational wave function (Gutzwiller, 1963 Gutzwiller, 1964 Gutzwiller, 1965). This 
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Fig. 1.9 Power-law distributions of quasiparticle weiglits Z in tiie lialf-filled disordered Hub- 
bard model for different interaction strengths {U/Uc (W) = 0.6, 0.8, 0.9, 0.92, 0.94, 0.97) show- 
ing characteritic power-law behavior P (Z) ~ Z°'~^: as (7 — > Uc, a — >■ and the curves become 
increasingly singular at small Z. The inset shows the weak dependence on the lattice size. 



From reference ( Andrade et al, 2009 a I 



is not too surprising, as the Kotliar-Ruckenstein theory, when apphed to the Mott 
transition, is known to be equivalent to the GutzwiUer wave function approach. This 
kind of approach is known to be able to capture the low-energy features, such as the 
disappearance of the quasiparticle peak as the transition is approached, unlike the HI 
approximation. In the impurity problem language, the low-energy sector is described 
by two parameters, like in the infinite-C/ case (see Section 1.3.2): the quasiparticle 
weight Z, well known from Fermi liquid theory, which determines the width of the 
Kondo resonance (essentially the Kondo temperature) and the resonance position e, 
which measures its shift from the chemical potential. It is important to notice that in 
the clean lattice, Z also determines the effective carrier mass {m/m* ~ Z), a feature 
unique to cases in which the self-energy only depends on the frequency. Indeed, as the 
interaction strength is tuned to its critical value U Uc, Z ^ Uc — U 0, signal- 
ing the transmutation of the itinerant carriers into localized magnetic moments. In 
the statDMFT description, these local quantities vary from site to site, Z ^ Zi and 
£ — >■ Ei, and their distributions were thoroughly analyzed. Surprisingly, their critical 
behaviors were found to be very dissimilar. 

As the transition is approached, which now happens at a disorder-dependent critical 
interaction U c (W^), a-U — ^ 0, just like in the conventional Gutzwiller-Brinkman-Rice 
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scenario. However, the quasiparticle weight distribution P (Z) becomes increasingly 
broader as ?7 — > J7c {W). In fact, the typical value Ztyp = exp((lnZ)) — 0, whereas 
the mean value remains finite, indicating that although almost all sites become local 



moments, some remain empty or doubly occupied (Aguiar et al., 2006). Besides, the 



Z distribution acquires a generic power-law shape, as in other Griffiths phases (see 

P(Z)-Z"~\ (1.60) 



Fig. (1.91) 



This is very similar to the Bethe lattice studies of Section 1.3.2 and like in those cases, 
these power laws generate a singular thermodynamic response, see Eq. (1.551. Never- 



theless, whereas before we had a disorder-driven Anderson-type transition, here this 
generic behavior is found in the proximity of the interaction-driven Mott transition 
for fixed disorder strength^ clearly showing the amplifying effects of electronic correla- 
tions. The similarity to the quantum Griffiths scenario of magnets is not fortuitous. 
The low-Z values which dominate the thermodynamics occur in exponentially rare 
regions of suppressed disorder. Finally, we note that just like in other known quantum 
Griffiths phases, this electronic Griffiths phase seems to be tied to a phase transition 
characterized by an infinite randomness fixed point: we find that, up to the numerical 
uncertainty, a — >■ as C/ — >■ ?7c {W). 

Interestingly, correlations have the opposite effect on the distribution of resonance 
positions e^. Indeed, the width of the e distribution decreases as C/ — >■ Uc{W) (see 



Fig. (l.lOl). This is easily understood from the pinning of the Kondo resonances to 



the chemical potential, as already noted within DMFT (Tanaskovic et al., 2003) (see 
Section 1.1.5). Just like in DMFT, this leads to a strong disorder screening effect, 



although, unlike in DMFT, here the screening effect is not perfect: a small amount of 
disorder seems to survive as the transition is approached. 

The disordered Hubbard model has also been investigated very recently in two di- 



mensions with statDMFT and iterative perturbation theory as impurity solver ( Semm 



ler et al., 2010a). The disorder model used, however, was not the usual uniform diago- 
nal disorder, but rather a combination of diagonal and off-diagonal disorder, together 
with disorder in the interaction term. This choice was intended to describe the so- 
called speckle disorder found in some set-ups of fermionic cold atoms loaded in optical 
lattices. The phase diagram was determined at half filling at zero as well as finite tem- 
peratures. A feature specific to this kind of disorder is the fact that it is unbounded 
and arbitrarily large values of on-site potentials occur. As a result, long (exponen- 
tial) tails arise flanking the Hubbard bands which contribute to filling up a possible 
interaction-induced Hubbard gap. The main consequence of the presence of these tails 
is a partial suppression of the Mott insulating phase, in favor of a disordered strongly 
correlated metal phase. As a result, the Mott insulator and Anderson(-Mott) insulator 
phases do not share a phase boundary and are separated by a metallic phase. This is 



in contrast to the phase diagram obtained within TMT, see Section 1.2.2 



A single impurity in a strongly correlated host. The smallest element in an analysis 
of the effects of inhomogeneities in a crystal is a single isolated point-like impurity in 
an otherwise homogeneous host. In the case of a weakly correlated host, the expected 
effect is the formation of a potential scattering center, which can be usually be quite 
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Fig. 1.10 Distributions of the ratio of tiie local Kondo resonance position Si to the quasipar- 
ticle weight Zi, Vi = ii/Zi in the half-filled disordered Hubbard model showing a decreasing 
effective disorder strength, as measured by Si, with increasing interactions. The inset shows 
the weak dependence on the lattice size (full line for L = 20 and dashed line for L = 50). 



From reference ( Andrade et al, 200961 



simply described at low energies by a set of phase shifts with a negligible temperature 
dependence. If the host is a weakly correlated metal, the localized impurity generates 
a radially oscillatory disturbance of the host charge. The oscillation is characterized by 
a wave vector determined by the extremal radii of the Fermi surface and an envelope 
which decays as the inverse of the distance to the impurity raised to the power d in d 
dimensions 

cos(2fcFr) 

5n{x) -. . (1.61) 



These are known as Friedel oscillations (Ruderman and Kittel, 1954 Kasuya, 1956 



Yosida, 19571. A dilute collection of these impurities can then be treated as indepen- 
dent and their effects on transport properties is readily computed by conventional 
techniques, such as the Boltzmann equation. 

The situation can be very different if the host is a strongly correlated system. The 
local inhomogeneity can disrupt the delicate balance of its host and lead to a spatially 
dependent pattern which has to be computed in a fully self-consistent manner. For 
example, local moments or antiferromagnetic ordering can be induced in the vicinity 
of the impurity (AUoul et al., 2009). Superconducting correlations, specially of the 



unconventional type (e. g., d-wave), also lead to characteristic spatial patterns in the 
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vicinity of impurities, which can be probed by scanning tunnehng microscopy (STM) 
techniques (Balatsky et al., 2006). Besides, strong correlations typically generate very 



low energy scales. Therefore, one should expect a non-negligible energy and tempera- 
ture dependence with possibly observable effects. Indeed, STM studies have revealed 
a non-trivial interplay between spatial fluctuations and the energy dependence of the 



local density of states in the cuprate superconductors (McElroy et al., 2005). In this 



case, the spectral function, probed by STM, is hardly affected by disorder up to the 
superconducting gap energy, whereas much stronger spatial fluctuations are observed 
at higher energies. 

In general, the localized imperfection can be viewed as a source probe coupled to 
all wave vectors of the host charge. Therefore, if this coupling can be treated within 
linear response theory, the relevant quantity to keep track of is the host wave-vector- 
dependent charge susceptibility. The vicinity to the Mott metal-insulator transition 
causes the charge susceptibility to be strongly suppressed as the system becomes in- 
creasingly more localized. Therefore, we would expect considerable changes in the 
spatial pattern of the Friedel oscillations in a strongly correlated metal. The stat- 
DMFT is the tool of choice to study this situation. A single non-magnetic impurity in 



a half-filled Hubbard model has been studied using just such a tool (Andrade et al, 



2010). A great deal of analytical insight can then be obtained through an expansion 



in the disorder strength to first non-trivial order, using the Kotliar-Ruckenstein slave- 



boson impurity solver (Kotliar and Ruckenstein, 1986 1. Note that the procedure is fully 



non-perturbative in the interaction strength. In the weak coupling limit, U <^ D, the 
approach has been shown to be fully equivalent to the Hartree-Fock approximation, 
which predicts a static self-energy 



Sj {iuJn) — Un 



(0) 



where n^-"'* is the non-interacting charge density 



(0) 1 



(1.62) 



(1.63) 



where n|^'is the real-space static Lindhard polarization function and V is the potential 
of the single impurity at site 0. This contains the usual Friedel oscillation pattern of 



Eq. (1.61) 



However, as the interaction strength is tuned to be close to the Mott localization 
limit Uc, a very different situation arises. In this case, the charge density pattern decays 



as (Andrade et al, 2010) 



no 



n(o) 



V (U^ Uc) . 



(1.64) 



Two features of this result stand out. First, the spatial pattern is determined by the 
inverse of the static Lindhard function. Surprisingly, however, up to logarithmic cor- 
rections, the inverse static Lindhard function behaves very much the same as the usual 
Lindhard function. Second, the amplitude of the strongly-correlated Friedel oscillations 
is significantly suppressed as U ^ Uc down by the third power of the parametric dis- 
tance to the transition. As a result, the charge perturbation dies out significantly as 
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Fig. 1.11 Friedel oscillation patterns of charge oscillations generated by a localized non- mag- 
netic impurity in a half-filled Hubbard model. The long-ranged oscillations of Srn = — 1 in 
the non-interacting case (top figure) are strongly suppressed in the strongly correlated regime 
(bottom figure, m/m* = 0.3). From reference (Andrade et al., 2010|. 



one recedes from the impurity, implying a short "healing length" (see Fig. 1.11 ). Com 



parison with numerical solutions of the statDMFT equations shows that the analytical 
expressions are quite accurate even for \V\ ^ D. 

The scattering T-matrix, which governs the transport properties, can also be com- 
puted with the same technique. We find its real space expression to be 
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T, = 



f/n 



(1.65) 



at small U and 



Uc-U 



n(o) 



iO 



(1.66) 



as one approaches the Mott transition. Here, the point to note is once again the strong 
suppression of scattering in the critical region, a situation by now already familiar 
( "disorder screening" ) . 

The temperature dependence of the scattering is also very special. Indeed, it is 
known that the resistivity in the ballistic regime in two dimensions exhibits an anoma- 
lous, non- Fermi liquid linear in temperature dependence (Zala et al., 2001 1, due to the 



coherent scattering off the Friedel oscillations created by the dilute impurities. Since 
we find these oscillations to be strongly suppressed by correlations it is important to 
study if the anomalous behavior survives. Besides, strong interactions also give rise to 
large inelastic scattering effects, and it is unclear if these are strong enough to mask 



the anomalous temperature dependence. The full analysis shows that (Andrade et al., 



2010) 



The interaction suppression of scattering (disorder screening), Eq. (1.661, acts 



to weaken both elastic and inelastic contributions from impurities, but does not 
destroy the non-Fermi liquid, anomalous behavior. 
• However, inelastic scattering limits considerably the temperature range in which 
the linear in T behavior of the resistivity is observed. These effects come both from 
the Friedel-oscillation regions, close to the impurities, as well as from the bulk, 
impurity-free regions of the system. It is found that, although both contributions 
are deleterious to the anomalous behavior, it is the bulk contribution which always 
dominates. In practice, even for quite modest mass renormalizations, the linear 
in T regime is probably already unobservable: for m/m* ~ 0.6, for example, it is 
restricted to T < IQ-^Tp, where Tp ~ D. 

Both effects above seem to make the observation of the anomalous resistivity in two 
dimensions very unlikely in the strongly correlated regime. 

In more general terms, the strategy of expansion in weak disorder strength outlined 
above, which is quite general, seems to be a good avenue for exploration of strongly 
correlated spatially inhomogeneous systems. 

1.3.4 Inhomogeneous systems with slab geometries 



Although the local approach outlined in Sections 1.3.1 1.3.2 and 1.3.3 has been de- 



scribed in the context of random systems, its microscopic formulation (Dobrosavljevic 



and Kotliar, 1997 Dobrosavljevic and Kotliar, 1998 ) is clearly applicable in any spa- 



tially non- uniform situation, whether random or not. Indeed, in the past several years 
this method has been applied to a number of non-random, yet spatially non-uniform 
systems. Because of this conceptual similarity, we chose to describe it here under the 
statDMFT "umbrella", although the systems have nothing "statistical" in them. A 
particular focus has been put on systems composed of different materials arranged 
in a slab geometry. Motivation for these investigations has come both from advances 
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in growth techniques of epitaxial heterostructures of complex materials, particularly 



transition metal oxides (Ohtomo et al., 2002), and also from the use of surface sensitive 



probes of strongly correlated materials, such as photoemission experiments ( Sckiyama 



et al., 2000 


Mo et al, 2003 


Rodolakis et al., 2009 


I. A number of different spatial 
ials have been considered: semi- 
iing planar boundary ("strongly 


arrangements of complex, strc 
infinite strongly correlated sys 
correlated surface physics") ( 


mgly correlated mater 
items with a free-stanc 


Potthoff and Nolting, 1999 Liebsch, 2003 Schwieger 


et al, 2003 


Borghi et al, 2009 


Borghi et al, 2010 


), planar interfaces between two 



semi-infinite materials (Helmes et al, 2008), heterostructures (Okamoto and Millis, 



2004a Okamoto and Millis, 2004&), semi-infinite metallic leads sandwiching barriers 



of strongly correlated metals or insulators (Freericks, 2004 Chen and Freericks, 2007 



Zenia et al, 2009 


Borghi et al, 2010 


). s 


also been performed ( 


Snoek et al, 2008) 



Studies of other types of inhomogeneities have 



The abrupt change of chemical properties at the planes of surfaces or interfaces 
can lead to novel phenomena, some of which have been identified and characterized 
by the local approach we are describing here. This a fast moving area of research and 
we will not attempt to cover the richness of behavior that has been and continues to 
be uncovered. We will confine ourselves to a few important findings which serve to 
convey the flavor of these phenomena. 

The interface between a metal and a Mott insulator. Metal-insulator and metal- 
semiconductor interfaces have been produced and studied for many years. Phenomena 
such as band bending due to charge rearrangement are well known. However, the 
question of what happens when a metal is grown epitaxially onto a Mott insulator has 
received much less attention. A first attempt to elucidate this question has been made 
by considering an ideal interface between particle-hole symmetric metallic and Mott 



insulating systems (Helmes et al, 2008 Borghi et al, 2009 Borghi et al, 2010). In 



the absence of antiferromagnetism and within the DMFT picture, the spins of a bulk 
Mott insulator are not quenched because there are no conduction electrons available at 
the Fermi level with which they can form a singlet. When a Mott insulator is brought 
into contact with a metal, however, the layers which are closest to the interface can 
have access to the metallic carriers (through tunneling) and be quenched by means 
of the Kondo effect. This has been dubbed the "Kondo proximity effect". Evidently, 
the metallization and Kondo quenching processes are established in a spatially smooth 
fashion as one enters the insulator, thus creating a strongly correlated metallic surface 
layer. The metallic character can be characterized by means of the by now familiar 
quasiparticle weight Z (x), which is here a spatially varying quantity that depends 
on the distance x from the interface. It is found that Z (x) decays exponentially as a 
function of x 



Z{x) - exp(-x/^) , 



(1.67) 



with a characteristic length scale ^ which sets the size of the metallic surface layer and 
is determined by the parametric distance to the Mott transition, ^ ~ (C/c ~ C^) ■ 
The 1/2 value of the critical exponent is expected from the mean- field character of the 
approach. At Mott criticality, this characteristic length diverges and 
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(1.68) 



The numerical value of the constant A above is found to be extremely small (^ 0.008), 
showing the Kondo induced penetration of the metal into the Mott insulator to be 
extremely ineffective. The small values of the quasiparticle residues, which govern the 
local Fermi energy scale, also lead to a strong dependence on temperature, energy 
or applied voltage. Very small values of these perturbations are able to completely 



destroy the metallic behavior, leading to the concept of a "fragile" Fermi liquid ( Zenia 



et al., 2009) 



Surface dead layer. Another important phenomenon occurs at the free planar sur- 



face bounding a strongly correlated metal (Borghi et al., 2009 Borghi et al., 2010). 
Here again the proximity to the interface with the vacuum induces variations of the 
quasiparticle residue Z. Interestingly, since the outermost layer electrons cannot tun- 
nel further out of the material, they have less of a chance to hybridize and stabilize the 
metallic behavior. As a consequence, there is a strong suppression of the quasiparticle 
weight close to the open surface. In effect, the outermost electrons are a very "fragile" 
Fermi liquid if not completely Mott localized, thus creating a surface layer of almost 
Mott insulating character, the so-called "dead layer". Once again the characteristic 

— 1 /2 

width of this layer also depends on the proximity to criticality as ^ ~ (C/ — Uc) ■ 

This has important consequences for the interpretation of photoemission experi- 
ments. For a long time, a quasiparticle peak had been sought in the photoemission 
spectra of strongly correlated materials with little success. While there have been many 
attempts to explain this mystery, it is clear now that one expects on general grounds 
that the quasiparticle peak width (which is proportional to Z) becomes extremely 
narrow close to the surface. Since conventional photoemission spectra only probe the 
outer surface layers of the material, it is not too surprising that it has been difficult 
to detect a well-formed quasiparticle peak in the past. More recently, however, higher- 
energy photons have been used, which are able to penetrate deeper into the compound 
and thus probe the behavior more typical of the bulk. As expected from the theoretical 
results outlined above, the quasiparticle becomes much better defined with increasing 



incident photon energy (Sekiyama et al., 2000 Mo et al., 2003 Rodolakis et al., 2009) 



1.4 Glassy behavior of correlated electrons 

We discussed so far the nature of strongly inhomogeneous metallic phases, result- 
ing from the interplay of strong correlations and disorder. Such 'electronic Griffiths 
phases", which can be viewed as precursors to the metal-insulator transition, reflect 
the formation of rare regions with anomalously slow dynamics. The resulting non- 
Fermi-liquid behavior is generically characterized by power law anomalies, with non- 
universal, rapidly varying exponents. In contrast, many experimental data, especially 
in Kondo alloys, seem to show reasonably weak anomalies, close to marginal Fermi 
hquid behavior (Stewart, 2001). 
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1.4.1 Instability of the electronic Griffiths phases to spin-glass ordering 

Physically, it is clear what is missing from the theory. Similarly as magnetic GrilRths 



phases (Miranda and Dobrosavljevic, 2005 Vojta, 2006), the electronic Griffiths phase 
is characterized ( [Miranda and Dobrosavljevic, 2001| Tanaskovic et al., 2004[ ) by a broad 



distribution P (T^f) ~ (Tjf)" of local energy scales (Kondo temperatures), with the 
exponent a ^ W^^ rapidly decreasing with disorder W. At any given temperature, 
the local moments with Tnii) < T remain unscreened. As disorder increases, the num- 



ber of such unscreened spins rapidly proliferates. Within the existing theory (Miranda 



et al., 1996{ [Miranda et al., 1997a[ [Miranda and Dobrosavljevic, 2001] [Tanaskovid 



et al., 2004) these unscreened spins act essentially as free local moments and provide a 
very large contribution to the thermodynamic response. In a more realistic description, 
however, even the Kondo-unscreened spins are not completely free, since the metal- 
lic host generates long-ranged Ruderman-Kittel-Kasuya-Yosida (RKKY) interactions 
even between relatively distant spins. 

In a disordered metal, impurity scattering introduces random phase fluctuations 
in the usual periodic oscillations of the RKKY interaction, which, however, retains its 



power law form (although its average value decays exponentially ( Jagannathan et al., 



1988 Narozhny et al, 2001 )). Hence, such an interaction acquires a random amplitude 



Jij of zero mean but finite variance (Jagannathan et al., 1988 Narozhny et al., 2001) 



(4iR))-j^a- (1-69) 

As a result, in a disordered metallic host, a given spin is effectively coupled with 
random but long range interactions to many other spins, often leading to spin-glass 
freezing at the lowest temperatures. How this effect is particularly important in Grif- 



fiths phases can also be seen from the mean-field stability criterion (Bray and Moore, 



1980 ) for spin glass ordering, which takes the form 

JxiociT) = 1. (1.70) 

Here, J is a characteristic interaction scale for the RKKY interactions, and xioc{T) is 
the disorder average of the local spin susceptibility. As we generally expect xioc{T) to 
diverge within a Griffiths phase, this arguments strongly suggests that in the presence 
of RKKY interactions such systems should have an inherent instability to finite (even 
if very low) temperature spin glass ordering. 

Similarly as other forms of magnetic order, the spin glass ordering is typically re- 
duced by quantum fiuctuations (e.g. the Kondo effect) which are enhanced by coupling 
of the local moments to itinerant electrons. Sufficiently strong quantum fluctuations 
can completely suppress spin-glass ordering even at T = 0, leading to a quantum crit- 
ical point separating a metallic spin glass from the conventional Fermi liquid ground 
state. As in other QCP's, one expects the precursors to magnetic ordering to emerge 
even before the transition is reached and produce non-Fermi liquid behavior within 
the corresponding quantum critical region. Since many systems where disorder-driven 
non-Fermi-liquid behavior is observed are not too far from incipient spin-glass order- 
ing, it is likely that these effects play an important role and should be theoretically 
examined in detail. 
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From the theoretical point of view, a number of recent works have examined the 
general role of quantum fluctuations in glassy systems and the associated quantum 
critical behavior. Most of the results obtained so far have concentrated on the behavior 
within the mean- field picture (i.e., in the limit of large coordination), where a consistent 
description of the quantum critical point (QCP) behavior has been obtained for several 



models. In a few cases (Read et al., 1995), corrections to mean-field theory have been 
examined, but the results appear inconclusive and controversial at this time. In the 
following, we briefly review the most important results obtained within the mean-field 
approaches. 

Quantum critical behavior in insulating and metallic spin glasses. 

• Ising spin glass in a transverse field 

The simplest framework to study the quantum critical behavior of spin glasses is 
provided by localized spin models such as the infinite-range Ising model in a transverse 
field (TFIM) with random exchange interactions Jij of zero mean and variance J"^ /N 
{N — > oo is the number of lattice sites) . 

Htfim = - E J'^j^t<^^ - r E (1-71) 

ij i 

In the classical limit (F = 0), this model reduces to the well-studied Sherrington- 



Kirkpatrick model (Mezard et al., 1986), where spins freeze with random orientations 
below a critical temperature Tsg(^ = 0) = J. Quantum fluctuations are introduced 
by turning on the transverse field, which induces up-down spin flips with tunneling 
rate '-^ F. As F grows, the critical temperature Tsg(^) decreases, until the quantum 
critical point is reached at F = Fj, « 0.731J, signaling the T = transition from a 



spin-glass to a quantum-disordered paramagnetic state (Fig. |1.12 ) 



Similarly as in DMFT theories for electronic systems, such infinite range models 
can be formally reduced to a self-consistent solution of an appropriate quantum im- 
purity problem, as first discussed in the context of quantum spin glasses by Bray and 



Moore (Bray and Moore, 1980). Early work quickly established the phase diagram 



( Dobrosavljevic and Stratt, 1987) of this model, but the dynamics near the quan- 
tum critical point proved more difficult to unravel, even when the critical point is 
approached form the quantum-disordered side. Here, the problem reduces to solving 
for the dynamics of a single Ising spin in a transverse field, described by an effective 



Hamiltonian (Miller and Huse, 1993) of the form 



H = J2 J J dTdT'a'{T)x{T - r')^T^(T') +T j dra^r). 

Physically, the interaction of the considered spin with the spin fluctuations of its 
environment generates the retarded interaction described by the "memory kernel" 
x(t — t'). An appropriate self-consistency condition relates the memory kernel to the 
disorder-averaged local dynamical susceptibility of the quantum spin 



x(T-T') = (ra-(T)a-(T')). 
A complete solution of the quantum critical behavior can be obtained, as first estab- 



lished in a pioneering work by Miller and Huse ( Miller and Huse, 1993 1 . These authors 
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behavior for spin glasses. The parameter r, which measures the quantum fluctuations, can 
represent the transverse field for localized spin models or the Fermi energy in metallic spin 
glasses. 

have set up a diagrammatic perturbation theory for the dynamic susceptibility, show- 
ing that the leading loop approximation already captures the exact quantum critical 
behavior, as the higher order corrections provide only quantitative renormalizations. 
The dynamical susceptibility takes the general form 

where the local static susceptibility Xo remains finite throughout the critical regime, 
and the spin excitations exist above a gap 

A - (r/|lnr|)^/^ 

which vanishes as the transition is approached from the paramagnetic side (here r = 
(r — Fc) /Fc measures the distance from the critical point). 

• The quantum spin-glass phase and the replicon mode 



The validity of this solution was confirmed by a generalization (Ye et ai, 1993) 
to the M-component rotor model (the Ising model belongs to the same universality 
class as the M — 2 rotor model), which can be solved in closed form in the large M 
limit. This result, which proves to be exact to all orders in the 1/M expansion, could 
be extended even to the spin glass phase, where a full replica symmetric solution was 
obtained. Most remarkably, the spin excitation spectrum remains gapless (A = 0) 
throughout the ordered phase. Such gapless excitations commonly occur in ordered 
states with broken continuous symmetry, but are generally not expected in classical 
or quantum models with a discrete symmetry of the order parameter. In glassy phases 
(at least within mean- field solutions) , however, gapless excitations generically arise for 
both classical and quantum models. Here, they reflect the marginal stability (Mezard 
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et al., 1986) found in the presence of replica symmetry breaking, a phenomenon which 



reflects the high degree of frustration in these systems. The role of the Goldstone mode 
in this case is played by the so-called "replicon" mode, which describes the collective 
low energy excitations characterizing the glassy state. 

A proper treatment of the low energy excitations in this regime requires special 
attention to the role of replica symmetry breaking (RSB) in the T limit. The 



original work (Ye et al., 1993) suggested that RSB is suppressed at T = 0, so that 



the simpler replica symmetric solution can be used at low temperatures. Later work 



(Georges et al, 2001), however, established that the full RSB solution must be consid- 
ered before taking the T — >■ limit, and only then can the correct form of the leading 
low temperature corrections (e.g., the linear T-dependence of the specific heat) be 
obtained. 

• Physical content of the mean-field solution 



In appropriate path-integral language (Dobrosavljevic and Stratt, 1987 Ye et al., 



1993), the problem can be shown to reduce to solving a one dimensional classical 



Ising model with long-range interactions, the form of which must be self-consistently 
determined. Such classical spin chains with long range interactions in general can be 



highly nontrivial. Some important examples are the Kondo problem (Anderson and 



Yuval, 1969 Yuval and Anderson, 1970 Anderson et al., 1970), and the dissipative 



two-level system (Leggett et al., 1987), both of which map to an Ising chain with 
1/r^ interactions. Quantum phase transitions in these problems correspond to the 



Kosterlitz-Thouless transition found in the Ising chain (Kosterlitz, 1976), the descrip 



tion of which required a sophisticated renormalization-group analysis. Why then is 
the solution of the quantum Ising spin glass model so simple? The answer was pro- 
vided in the paper by Ye, Sachdev, and Read (Ye et al, 19931, which emphasized 



that the critical state (and the RSB spin-glass state) does not correspond to the crit- 
ical point, but rather to the high-temperature phase of the equivalent Ising chain, 
where a perturbative solution is sufficient. In Kondo language, this state corresponds 
to the Fermi-liquid solution characterized by a finite Kondo scale, as demonstrated 



by a quantum Monte-Carlo calculation of Rozenberg and Grempel (Rozenberg and 



Grempel, 1998), which also confirmed other predictions of the analytical theory. 



From a more general perspective, the possibility of obtaining a simple analytical 
solution for quantum critical dynamics has a simple origin. It follows from the fact 
that all corrections to Gaussian (i.e. Landau) theory are irrelevant above the upper 



critical dimension, as first established by the Hertz-Millis theory (Hertz, 1976 Millis 



1993 ) for conventional quantum criticality. The mean-field models become exact in 



the limit of infinite dimensions, hence the Gaussian solution of Refs. (Miller and Huse, 



1993 Ye et al., 1993 ) becomes exact. Leading corrections to mean-field theory for rotor 



models were examined by an £-expansion below the upper critical dimension dc = 8 
for the rotor models by Read, Sachdev, and Ye, but these studies found run-away 
fiows, presumably indicating non-perturbative effects that require more sophisticated 
theoretical tools. Most likely these include Griffiths phase phenomena controlled by 



the infinite randomness fixed point, as already discussed in Section 1.3 



• Metallic spin glasses 
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A particularly interesting role of the low-lying excitations associated with the spin- 
glass phase is found in metallic spin glasses. Here the quantum fluctuations are pro- 
vided by the Kondo coupling between the conduction electrons and local moments, 
and therefore can be tuned by controlling the Fermi energy in the system. The situa- 
tion is again the simplest for Ising spins where an itinerant version of the rotor model 
of Sengupta and Georges ( Sengupta and Georges, 1995 1 can be considered, and similar 



results have obtained for the "spin-density glass" model of Sachdev, Read, and Opper- 



mann (Sachdev et ai, 1995 1. The essential new feature in these models is the presence 



of itinerant electrons which, as in the Hertz-Millis approach (Hertz, 1976 Millis, 1993), 



have to be formally integrated out before an effective order-parameter theory can be 
obtained. This is justified provided that the quasiparticles remain well defined at the 
quantum critical point, i.e. that the quasiparticle weight Z ~ Tr- remains finite and 
the Kondo effect remains operative. The validity of these assumptions is by no means 
obvious, and led to considerable controversy before a detailed quantum Monte Carlo 



solution of the model became available (Rozenberg and Grempel, 1999), confirming 
the proposed scenario. 

Under these assumptions, the theory can again be solved in closed form, and we 
only quote the principal results. Physically, the essential modification is that the pres- 
ence of itinerant electrons now induces Landau damping, which creates dissipation 
for the collective mode. As a result, the dynamics is modified, and the local dynamic 
susceptibility now takes the following form 



(1.72) 



The dynamics is characterized by the crossover scale uj* ^ r [r measures the distance to 
the transition) which defines a crossover temperature T* ^ uj* separating the Fermi 
liquid regime (at T <C T*) from the quantum critical regime (at T ^ T*). At the 

1/2 

critical point xi'^n) — Xo + \^n\ , leading to non-Fermi liquid behavior of all physical 
quantities, which acquire a leading low-temperature correction of the T^^^ form. This 
is a rather mild violation of Fermi liquid theory, since both the static spin susceptibility 
and the specific heat coefficient remain finite at the QCP. A more interesting feature, 
which is specific to glassy systems, is the persistence of such quantum critical non- 
Fermi-liquid behavior throughout the metallic glass phase, reflecting the role of the 
replicon mode. 

Spin-liquid behavior, destruction of the Kondo effect by bosonic dissipation, and frac- 
tionaUzation. 

• Quantum Heisenberg spin glass and the spin-liquid solution 
Quantum spin glass behavior proves to be much more interesting in the case of 



Heisenberg spins, where the Berry phase term (Fradkin, 1991 ) plays a highly nontrivial 



role, completely changing the dynamics even within the paramagnetic phase. While 
the existence of a finite temperature spin-glass transition was established even in early 



work (Bray and Moore, 1980), solving for the details of the dynamics proved difficult 



until the remarkable work of Sachdev and Ye (Sachdev and Ye, 19931. By a clever use 



of large- iV methods, these authors identified a striking spin-liquid solution within the 
paramagnetic phase. In contrast to the nonsingular behavior of Ising or rotor quantum 
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spin glasses, the dynamical susceptibility now displays a logarithmic singularity at low 
frequency. On the real axis it takes the form 

X{i^) ln(l/ \uj\) + i^sgn{uj). 

A notable feature of this solution is that it is precisely of the form postulated for 



"marginal Fermi liquid" phenomenology (Varma et al., 1989) of doped cuprates. The 
specific heat is also found to assume a singular form C ^ VT, which was shown 



(Georges et al., 2001) to reflect a nonzero extensive entropy if the spin liquid solution 
is extrapolated to T = 0. Of course, the spin-liquid solution becomes unstable at a 
finite ordering temperature, and the broken symmetry state has to be examined to 
discuss the low temperature properties of the model. 

Subsequent work ( Georges et al, 2001| ) demonstrated that this mean-field solution 
remains valid for all finite N and generalized the solution to the spin-glass (ordered) 
phase. A closed set of equations describing the low temperature thermodynamics in the 
spin glass phase was obtained, which was very recently re-examined in detail ( [Camjayi 



and Rozenberg, 2003 1 to reveal fairly complicated behavior. 



• Metallic Heisenberg spin glasses and fractionalization 

Even more interesting is the fate of this spin liquid solution in itinerant systems, 
where an additional Kondo coupling is added between the local moments and the 
conduction electrons. The mean-field approach can be extended to this interesting 



situation by examining a Kondo-Heisenberg spin glass model (Burdin et al., 2002) 
with the Hamiltonian 

Hkh = -tJ2 (^-S- + H- c.) + 51 • s, + ^ J,,S, • S,. (1.73) 

{ij)cr i (ij) 

1 /2 

In the regime where the scale of the RKKY interaction J — (Jfj) is small 
compared to the Kondo coupling Jk , one expects Kondo screening to result in standard 
Fermi liquid behavior. In the opposite limit, however, the spin fluctuations associated 
with the retarded RKKY interactions may be able to adversely affect the Kondo 
screening, and novel metallic behavior could emerge. This intriguing possibility can 
be precisely investigated in the mean-field (infinite range) limit, where the problem 
reduces to a single- impurity action of the form ( [Burdin et al., 20"02 ) 



„ JO 



^ dr r dr'cUr)G,{r-T'Kir') 



Jk / dTS(r) • s(t) 
Jq 

/ dT / dT'x(r-T')S(r).S(T'). (1.74) 
JO Jo 



j2 rP 
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RKKY coupling J 



Fig. 1.13 Phase diagram of the Bose-Fermi Kondo model in the presence of a sub-Ohmic 
bosonic bath ( |Si and Smith, 1996| [Sengupta, 2000| |Zhu and Si, 2002| [Zarand and Demler,| 
2002 I . Kondo screening is destroyed for sufficiently large dissipation (RKKY coupling to spin 
fluctuations). 



Such a single- impurity action (1.74 1 describes the so-called Bose-Fermi Kondo 



(BFK) impurity model (Si and Smith, 1996 Sengupta, 2000 Zhu and Si, 2002 Zarand 



and Demler, 2002 1 where, in addition to the coupling to the fermionic bath of conduc- 



tion electrons, the Kondo spin also interacts with a bosonic bath of spin fluctuations, 
with local spectral density x i^n)- Because the same BFK model also appears in "ex- 
tended" DMFT theories ( Smith and Si, 2000 1 of quantum criticality in clean systems 



(Si and Smith, 1996 Sengupta, 2000 Si et al., 2001 Zhu and Si, 2002 Zarand and 



Demler, 2002), its properties have been studied in detail and are by now well under- 
stood. 

In the absence of the RKKY coupling (J = 0), the ground state of the impurity 
is a Kondo singlet for any value of Jk 0. By contrast, when J > 0, the dissipation 
induced by the bosonic bath tends to destabilize the Kondo effect. For a bosonic bath 
of "Ohniic" form (x (w^) = Xo — C|a;„|), this effect only leads to a finite decrease 
of the Kondo temperature, but the Fermi liquid behavior persists. In contrast, for 
"sub-Ohmic" dissipation (x(i^n) = Xo ~ C \^n\^ *^ with e > 0) two different phases 
exist, and for sufficiently large RKKY coupling the Kondo effect is destroyed. The two 



regimes are separated by a quantum phase transition (see Fig. 1.13). 

Of course, in the considered Kondo lattice model with additional RKKY inter- 
actions, the form of the bosonic bath x i'^n) is self-consistently determined and can 
take different forms as the RKKY coupling J is increased. The model was analytically 
solved within a large N approach by Burdin et al. ( Burdin et al., 2002) , who calculated 



the evolution of the Fermi liquid coherence scale T* and the corresponding quasipar- 
ticle weight Z in the presence of RKKY interactions. Within the paramagnetic phase 
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both T*{J) and Z{J) are found to decrease with J until the Kondo effect (and thus 
the Fermi hquid) is destroyed at J = Jc « lOT^ (here is the J = Kondo tem- 
perature), where both scales vanish (Burdin et ai, 2002 Tanaskovic et ai, 2005). At 



T > T*{J) (and of course at any temperature for J > Jc) the spins effectively decou- 
ple from conduction electrons and spin liquid behavior, essentially identical to that of 
the insulating model, is established. Thus, sufficiently strong and frustrating RKKY 
interactions are able to suppress Fermi liquid behavior, and marginal Fermi liquid be- 
havior emerges in a metallic system. This phenomenon, corresponding to spin-charge 
separation resulting from the destruction of the Kondo effect, is sometimes called 



"fractionalization" ( 


Coleman and Andrei, 1987 


2002 


Senthil et al, 2003 


Senthil et al., 2004 


1. 



Kagan et ai, 1992 Demler et al. 



cated as an appealing scenario for exotic phases of strongly correlated electrons, but 
with the exception of the described model, there are very few well established results 
and model calculations to support its validity. Finally, we should mention related work 



(Parcollet and Georges, 1999) on doped Mott insulators with random exchanges, with 
many similarities with the above picture. 

We should note, however, that this exotic solution is valid only within the param- 
agnetic phase, which is generally expected to become unstable to magnetic (spin glass) 
ordering at sufficiently low temperatures. Since fractionalization emerges only for suf- 
ficiently large RKKY coupling (in the large N model Jc ~ lOT^), while in general one 
expects magnetic ordering to take place already at J ~ (according to the famous 



Doniach criterion (Doniach, 19771), one expects (Burdin et ai, 2002) the system to 
magnetically order much before the Kondo temperature vanishes. If this is true, then 
one expects the quantum critical behavior to be very similar to metallic Ising spin 
glasses, i.e., to assume the conventional Hertz- Millis form, at least for the mean- field 
spin glass models we discussed here. The precise relevance of this paramagnetic spin 
liquid solution thus remains unclear, at least for systems with weak or no disorder in 
the conduction band. 

• Fractionalized two-fluid behavior of electronic Griffiths phases 

The situation seems more promising in the presence of sufficient amounts of dis- 
order, where the electronic Griffiths phase forms. Here the disorder generates a very 
broad distribution of local Kondo temperatures, making the system much more sen- 
sitive to RKKY interactions. This mechanism has recently been studied within an 



extended DMFT approach (Tanaskovic et al., 2005), which is able to incorporate both 



the formation of the Griffiths phase and the effects of frustrating RKKY interactions 
leading to spin-glass dynamics. At the local impurity level, the problem is still reduced 
to the Bose-Fermi Kondo model, but the presence of conduction electron disorder qual- 
itatively modifies the self-consistency conditions determining the form of x (wn)- 

To obtain a sufficient condition for decoupling, we examine the stability of the 
Fermi liquid solution, by considering the limit of infinitesimal RKKY interactions. To 
leading order we replace 



X(t^„) > Xo (t^n) = Xi^^n, J = 0), 



and the calculation reduces to the "bare model" of Ref. (Tanaskovic et al., 2005) 
this case, P (Tk) ~ T^"\ where a ~ and 



In 
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X0K)~ J dTKPiTK)xi^n,TK)-Xo{0)~Co\LOn\^~\ (1.75) 

where e = 2 — a. Thus, for sufficiently strong disorder (i.e. within the electronic 
Griffiths phase), even the "bare" bosonic bath is sufficiently singular to generate de- 
coupling. The critical value of W will be modified by self-consistency, but it is clear 
that decoupling will occur for sufficiently large disorder. 

Once decoupling is present, the system is best viewed as composed of two fluids, 
one made up of a fraction n of decoupled spins, and the other of a fraction (1 — n) 
of Kondo screened spins. The self-consistent xi^n) acquires contributions from both 
fluids 

X (wn) = nxdc (wn) + {l-n) Xs (w„) . (1.76) 



A careful analysis (Tanaskovic et al., 2005) shows that, for a bath characterized by an 
exponent e 

Xdc {i^n) - Xdc (0) - C ; (1.77) 
Xs (oJn) - Xs (0) - C |^„|i-(2-^-^/'') , (1.78) 

where v = ly (e) is a critical exponent governing how the Kondo scale vanishes at the 
quantum critical point of the Bose-Fermi model. Since v > 0, the contribution of the 
decoupled fluid is more singular and dominates at lower frequencies. Self-consistency 



then yields e = 1, as in the familiar spin liquid state of Sachdev and Ye (Sachdev 



and Ye, 1993). For e — 1, the local susceptibility is logarithmically divergent (both in 



ijjn and T). This does not necessarily mean that the bulk susceptibility, which is the 



experimentally relevant quantity, behaves in the same manner (ParcoUet and Georges, 



1999). More work remains to be done to determine the precise low temperature form 
of this and other physical quantities and to assert the relevance of this mechanism for 
specific materials. 

As in the case where conduction electron disorder is absent, the spin liquid state 
is unstable towards spin-glass ordering at sufficiently low temperatures. However, nu- 



merical estimates for the Griffiths phase model (Tanaskovic et al., 2005) suggest a sur- 
prisingly wide temperature window where the marginal behavior should persist above 
the ordering temperature. Fig. |1.14| represents the predicted phase diagram of this 
model. For weak disorder the system is in the Fermi liquid phase, while for W > Wc 
the marginal Fermi liquid phase emerges. The crossover temperature (dashed line) 
delimiting this regime can be estimated from the frequency up to which the logarith- 
mic behavior of the local dynamical susceptibility xi^^) is observed. The spin glass 



phase, obtained from Eq. (1.70), appears only at the lowest temperatures, well below 



the marginal Fermi liquid boundary. Interestingly, recent experiments (MacLaughlin 



et al., 2001) have indeed found evidence of dynamical spin freezing in the milliKelvin 
temperature range for the same Kondo alloys that display normal phase non-Fermi- 
liquid behavior in a much broader temperature window. 

The two-fluid phenomenology of the disordered Kondo lattice we have described 
above is very reminiscent of earlier work on the clean Kondo lattice, where the conduc- 
tion electrons effectively decouple from the local moments, the latter forming a spin 
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Fig. 1.14 Phase diagram of tiie electronic Griffiths phase model with RKKY interactions 
( [Tanaskovic et al., 2005[ ) . The inset shows the fraction of decoupled spins as a function of the 
disorder strength W . 



liquid state 


(Coleman and Andrei, 1987 


Kagan et al, 1992 


Demler et al, 2002 


Senthil 


et al, 2003 


Senthil et al, 2004 


). The major difference between the results presented in 



this Section and these other cases is that here local spatial disorder fluctuations lead to 
an inhomogeneous coexistence of the two fluids, as each site decouples or not from the 
conduction electrons depending on its local properties. The discussed mean-field mod- 
els should be considered as merely the first examples of this fascinating physics. The 
specific features of the spin liquid behavior that was obtained from these models may 
very well prove to be too restrictive and perhaps even inaccurate. For example, the 
specific heat enhancements may well be overestimated, reflecting the residual T = 
entropy of the mean-field models. Nevertheless, the physics of Kondo screening being 
destroyed by the interplay of disorder and RKKY interactions will almost certainly 
play a central role in determining the properties of many non-Fermi-liquid systems, 
and clearly needs to be investigated in more detail in the future. 

1.4.2 Electron glass 

Another aspect of disordered interacting electrons poses a fundamental problem. Very 
generally. Coulomb repulsion favors a uniform electronic density, while disorder fa- 
vors local density fluctuations. When these two effects are comparable in magnitude, 
one can expect many different low energy electronic configurations, i.e. the emergence 
of many metastable states. Similarly as in other "frustrated" systems with disorder, 
such as spin glasses, these processes can be expected to lead to glassy behavior of the 
electrons, and the associated anomalously slow relaxational dynamics. Indeed, both 



theoretical ( 


Davies et al, 1982 


|Pollak and Hunt, 1991 


) and experimental ( 


Ben-Chorin 


et al, 1993||Ovadyahu and Po 


lak, 19971 


Bogdanovich and Popovic, 2002| 


Jaroszyhski 


et al, 2002 


1 work has found evidence of such behavior deep on the insulating side 



of the transition. However, at present very little is known as to the precise role of 
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such processes in the critical region. Nevertheless, it is plausible that the glassy freez- 
ing of the electrons must be important, since the associated slow relaxation clearly 
will reduce the mobility of the electrons. From this point of view, the glassy freezing 
of electrons may be considered, in addition to the Anderson and the Mott mecha- 
nism, as a third fundamental process associated with electron localization. Interest in 
understanding the glassy aspects of electron dynamics has experienced a genuine re- 
naissance in the last few years, primarily due to experimental advances. Emergence of 
many metastable states, slow relaxation and incoherent transport have been observed 
in a number of strongly correlated electronic systems. These included transition metal 
oxides such as high Tc materials, manganites, and ruthenates. Similar features have 
recently been reported in two-dimensional electron gases and even three-dimensional 
doped semiconductors such as Si:P. 

• Infinite-dimensional model of an electron glass 

The interplay of the electron-electron interactions and disorder is particularly ev- 
ident deep on the insulating side of the metal-insulator transition (MIT). Here, both 
experimental (Massey and Lee, 1996) and theoretical studies (Efros and Shklovskii, 



1975) have demonstrated that they can lead to the formation of a soft "Coulomb 



gap" , a phenomenon that is believed to be related to the glassy behavior ( Ben-Chorin 



et al, 1993 [Ovadyahu and PoUak, 1997| Bogdanovich and Popovic, 2002[ Jaroszynski 



et al., 2002) of the electrons. Such glassy freezing has long been suspected (Belitz and 



Kirkpatrick, 19951 to be of importance, but more recent theoretical work (DobrosavL 



jevic et al., 1997 



Chakravarty et al, 1999) has suggested that it may even dominate 



the MIT behavior in certain low carrier density systems. The classic work of Efros 



and Shklovskii (Efros and Shklovskii, 1975) has clarified some basic aspects of this 



behavior, but a number of key questions have remain unanswered. 



As a simplest example (Pastor and Dobrosavljevic, 1999 ) displaying glassy behavior 
of electrons, we focus on a simple lattice model of spinless electrons with nearest 
neighbor repulsion V in presence of random site energies and inter-site hopping t, 
as given by the Hamiltonian 



H 



<ij> 



-ttj + ei5ij)clcj + V 



<ij> 



(1.79) 



This model can be solved in a properly defined limit of large coordination number 
( Georges et a/., 1996[ ), where an extended dynamical mean-field theory (EDMFT) 



formulation becomes exact. We concentrate on the situation where the disorder (or 
more generally frustration) is large enough to suppress any uniform ordering. We then 
rescale both the hopping elements and the interaction amplitudes as tij — tij/^/z] 
Vij — >■ VijI yfz. As we will see shortly, the required fluctuations then survive even in 
the z oo limit, allowing for the existence of the glassy phase. Within this model: 



The universal form of the Coulomb gap (Efros and Shklovskii, 1975) proves to be a 

direct consequence of glassy freezing. 
The glass phase is identified through the emergence of an extensive number of 

metastable states, which in our formulation is manifested as a replica symme- 



try breaking instability ( Mezard et al., 1986, ) 
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* As a consequence of this ergodicity breaking (Mezard et ai, 1986), the zero- field 



cooled compressibility is found to vanish at T=0, suggesting the absence of screen- 



ing (Efros and Shklovskii, 1975) in disordered insulators 



* The quantum fluctuations can melt this glass even at T = 0, but the relevant 
energy scale is set by the electronic mobility and is therefore a nontrivial function 
of disorder. 

We should stress that although this model allows us to examine the interplay of 
glassy ordering and quantum fluctuations due to itinerant electrons, it is too simple 
to describe the effects of Anderson localization. These effects require extensions to 
lattices with finite coordination and will be discussed in the next section. 

For simplicity, we focus on a Bethe lattice at half filling and examine the z ^ oo 
limit. This strategy automatically introduces the correct order parameters and after 



standard manipulations (Dobrosavljevic and Kotliar, 1994) the problem reduces to a 



self-consistently defined single site problem, 
the form 



as defined by an the effective action of 



+ ^^'KMx(T,T')<5<(r')] + \v^Y.[ ly^'^^' <5nJ(r').(1.80) 



Here, we have used functional integration over replicated Grassmann fields ( DobrosavL 



jevic and Kotliar, 1994) c°(r) that represent electrons on site i and replica index a, 
and the random site energies are distributed according to a given probability dis- 
tribution P{si). The operators (5n"(T) = [c^"" {t)c^{t) — 1/2) represent the density 
fluctuations from half filling. The order parameters G{t~t'), x(t — r') and qab satisfy 
the following set of self-consistency conditions 



G(r -r')=j de.P{e,) < cfMc^T') 
X(r-r') = J de,P{e,) < 6nl''{T)Snt{r') >e//, 



Order parameters 



(1.81) 
(1.82) 
(1.83) 



In these equations, the averages are taken with respect to the effective action of 
Eq. (1.80). Physically, the "hybridization function" t^G{T — r') represents the single- 



particle electronic spectrum of the environment, as seen by an electron on site i. In 



particular, its imaginary part at zero frequency can be interpreted (Dobrosavljevic 



and Kotliar, 1994) as the inverse lifetime of the local electron and as such remains 



finite as long as the system is metallic. We recall (Dobrosavljevic and Kotliar, 1994) 
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that ioi V = these equations reduce to the famiUar CPA description of disordered 
electrons, which is exact for z — oo. The second quantity x(''" ~ ''"') represents an 
(interaction-induced) mode- coupling term that reflects the retarded response of the 
density fluctuations of the environment. Note that very similar objects appear in the 



well-known mode-coupling theories of the glass transition in dense liquids ( Cummins 



et ai, 19941. Finally the quantity qat (a 7^ b) is nothing but the familiar Edwards- 
Anderson order parameter qea- Its nonzero value indicates that the time averaged 
electronic density is spatially non-uniform. 

• Equivalent Infinite Range model 

From a technical point of view, a RSB analysis is typically carried out by focus- 
ing on a free energy expressed as a functional of the order parameters. In our Bethe 
lattice approach, one directly obtains the self-consistency conditions from appropriate 



recursion relations ( Dobrosavljevic and Kotliar, 1994) , without invoking a free energy 



functional. However, we have found it useful to map our z — oo model to another 
infinite range model, which has exactly the same set of order parameters and self- 
consistency conditions, but for which an appropriate free energy functional can easily 



be determined. The relevant model is still given Eq. (1.79), but this time with ran- 
dom hopping elements tij and random nearest- neighbor interaction Vij^ having zero 
mean and variances t"^, and V'^, respectively. For this model, standard manipulations 



(Dobrosavljevic and Kotliar, 1994) result in the following free energy functional 



In 



j de,P{e,) J Dcl''Dc'^ exp{~Seffii)} 



(1.84) 



with Seff{i) given by Eq. (1.80). The self-consistency conditions, Eqs. ( 1.81||1.83 ) then 
follow from 



= 6F/SGiT,T'y, = SF/SxiT,T'); - 6F/Sq,b 



(1.85) 



We stress that Eqs. ( 1.81|1.83 ) have been derived for the model with uniform 
hopping elements tij and interaction amplitudes Vij, in the z — oo limit, but the 
same equations hold for an infinite range model where these parameters are random 
variables. 

• The glass transition 

In our electronic model, the random site energies e,; play a role of static random 
fields. As a result, in the presence of disorder, the Edwards-Anderson parameter qea 
remains nonzero for any temperature, and thus cannot serve as an order parameter. 
To identify the glass transition, we search for a replica symmetry breaking (RSB) 
instability, following standard methods (jde Almeida and Thouless, 1978 Mezard and 
Young, 1992). We define Sqab — Qab — Q, and expand the free energy functional of Eq. 
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( |1.84 1 around the replica symmetric solution. The resulting quadratic form (Hessian 
matrix) has the matrix elements given by 



dqabdq. 



cd 



= {l3V)'^5acSbd + [< Sna{Ti)Snb{T2) >RS< Snc{T3)Snd{T4) >RS 



(1.86) 



r0 rP rP rP 

I / / / dTidT2dT3dT4[< 5na{Ti)dnb{T2)6nc(T3)5nd{T4) >rs], 
Jo Jo Jo Jo 

where the expectation values are calculated in the RS solution. Using standard manip- 



ulations (de Almeida and Thouless, 1978), and after lengthy algebra, we finally arrive 



at the desired RSB stability criterion that takes the form 



1 = 



iXloc{£z)] 



(1.87) 



Here, [...]dis indicates the a,vera,g6 over disorder a-nd \ioc 

(si) is the local compressibility, 

that can be expressed as 



XlociSi) 



d_l 



dr < Sni{T) >, 



(1. 



and which is evaluated by carrying out quantum averages for a fixed realization of 
disorder. The relevant expectation values have to be carried with respect to the full 
local effective action Seff{i) of Eq. (2), evaluated in the RS theory. In general, the 
required computations cannot be carried out in closed form, primarily due to the 
unknown "memory kernel" xi''' ~ ''"')■ However, as we will see, the algebra simplifies 
in several limits, where explicit expressions can be obtained. 

• Classical electron glass 

In the classical {t = 0) limit, the problem can be easily solved in closed form. We 
first focus on the replica symmetric (RS) solution and set Qab = q for all replica pairs. 
The corresponding equation reads 



+ 00 



dx 



^-'/^tanh^ 



-x{{f3Vfq+(j3W) 



(1. 



where we have considered a Gaussian distribution of random site energies of variance 
W"^. Note that the interactions introduce an effective, enhanced disorder strength 



Weff = + V^q, 



(1.90) 



since the frozen-in density fluctuations introduce an added component to the random 
potential seen by the electrons. As expected, g 7^ for any temperature when W j^O. 
If the interaction strength is appreciable as compared to disorder, we thus expect 
the resistivity to display an appreciable increase at low temperatures. We emphasize 
that this mechanism is different from Anderson localization, which is going to be 
discussed in the next section, but which also gives rise to a resistivity increase at low 
temperatures. 
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Next, we examine the instability to glassy ordering. In the classical {t = 0) limit 
Eq. (11. 871) reduces to 



+ 00 



dx 



-'/^cosh- 



(1.91) 



with Weff{q) given by Eq.(1.90). The resulting RSB instability line separates a low 
temperature glassy phase from a high temperature "bad metal" phase. At large disor- 
der, these expressions simplify and we find 



Tg 



1 V' 



W oo. 



(1.92) 



6V2^ W 

We conclude that Tq decreases at large disorder. This is to be expected, since in this 
limit the electrons drop into the lowest potential minima of the random potential. 
This defines a unique ground state, suppressing the frustration associated with the 
glassy ordering and thus reducing the glassy phase. It is important to note that for 
the well known de Almeida- Thouless (AT) line Tusb decreases exponentially in the 
strong field limit. In contrast, we find that in our case, Tq ~ 1/W decreases only slowly 
in the strong disorder limit. This is important, since the glassy phase is expected to 
be most relevant for disorder strengths sufficient to suppress uniform ordering. At the 
same time, glassy behavior will only be observable if the associated glass transition 
temperature remains appreciable. 

• The glassy phase 

To understand this behavior, we investigate the structure of the low-temperature 
glass phase. Consider the single-particle density of states at T=0, which in the classical 
limit can be expressed as 



P(£,i = 0) = ^^<5( 



(1.93) 



where = £i + ^ij^j ^-''^ the renormalized site energies. In the thermodynamic 
limit, this quantity is nothing but the probability distribution Pji{ef'). It is analogous 
to the "local field distribution" in the spin-glass models and can be easily shown to 
reduce to a simple Gaussian distribution in the RS theory, establishing the absence 
of any gap for T > Tq. Obtaining explicit results from a replica calculation in the 
glass phase is more difficult, but useful insight can be achieved by using standard 



simulation methods (Palmer and Pond, 1979 Pazmandi et ai, 1999) on our equivalent 



infinite-range model; some typical results are shown in Fig. 1.15 We find that as a 
result of glassy freezing, a pseudo-gap emerges in the single-particle density of states. 



reminiscent of the Coulomb gap of Efros and Shklovskii (ES) (Efros and Shklovskii, 
1975). The low energy form of this gap appears universal, 



p{e) K Ce^/V^, foi s < Eg; C = a = 1, 



(1.94) 



independent of the disorder strength again in striking analogy with the predictions 
of ES. To establish this result, we have used stability arguments very similar to those 
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Fig. 1.15 Single particle density of states in the classical {t — 0) limit at T = 0, as a 
function of disorder strength. Results are shown from a simulation on A'' = 200 site system, 
for W/V — 0.5 (thin line) and W/V = 1.0 (full line). Note that the low energy form of the 
gap takes a universal form, independent of the disorder strength W. The dashed line follows 
Eq. fTgii. 



developed for spin-glass (SG) models ( Pazmandi et al., 1999| ), demonstrating that 
the form of Eq. (1.94) represents an exact upper bound for p{e). For infinite- ranged 
SG models, as in our case, this bound appears to be saturated, leading to universal 
behavior. Such universality is often associated with a critical, self-organized state of 
the system. Recent work (Pazmandi et ai, 1999) finds strong numerical evidence of 
such criticality for SG models; we believe that the universal gap form in our case has 
the same origin. Furthermore, assuming that the universal form of Eq. (1.94) is obeyed 
immediately allows for an estimate of Tg{W). Using Eq. (1.87) to est imate the gap 
size for large disorder gives Tq ^ Eg ^ V^/W, in agreement with Eq. (1.92). 

The ergodicity breaking associated with the glassy freezing has important conse- 
quences for our model. Again, using the close similarity of our classical infinite range 
model to standard SG models (Mezard et al, 1986), it is not difficult to see that the 
zero-field cooled (ZFC) compressibility vanishes at T = 0, in contrast to the field-cooled 
one, which remains finite. Essentially, if the chemical potential is modified after the 
system is cooled to T = 0, the system immediately falls out of equilibrium and displays 
hysteretic behavior (Pazmandi et al., 1999) with vanishing typical compressibility. If 
this behavior persists in finite dimensions and for more realistic Coulomb interactions, 
it could explain the absence of screening in disordered insulators. 
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• Arbitrary lattices and finite coordination: mean-field glassy phase of the random- 
field Ising model 



The simplest theories of glassy freezing (Mezard et al., 1986) are obtained by exam 



ining models with random inter-site interactions. In the case of disordered electronic 
systems, the interactions are not random, but glassiness still emerges due to frustration 
introduced by the competition of the interactions and disorder. As we have seen for 



the Bethe lattice (Pastor and Dobrosavljevic, 1999), random interactions are gener 



ated by renormalization effects, so that standard DMFT approaches can still be used. 
However, one would like to develop systematic approaches for arbitrary lattices and in 
finite coordination. These issues already appear on the classical level, where our model 



reduces to the random-field Ising model (RFIM) (Nattermann, 1997). To investigate 
the glassy behavior of the RFIM (Pastor et al, 2002), a systematic approach that can 



incorporate short-range fluctuation corrections to the standard Bragg-Williams theory 
is the method of Plefka ( [Plefka, 1982| and Georges et al. ( [Georges et al, 1990[ ). This 
work has shown that: (i) corrections to even the lowest nontrivial order immediately 
result in the appearance of a glassy phase for sufficiently strong randomness; (ii) this 
low-order treatment is sufficient in the joint limit of large coordination and strong 
disorder; (iii) the structure of the resulting glassy phase is characterized by universal 
hysteresis and avalanche behavior emerging from the self-organized criticality of the 
ordered state. 



• Long-range Coulomb interactions and the Efros-Shklovskii gap 

So far we focused on the limit of large coordination, where the EDMFT approxima- 
tion is essentially exact and the resulting electron glass phase shows many similarities 



to the standard Parisi theory of spin glasses (Mezard et al., 1986). While it produced 



many appealing results, this mean-field approach has remained very controversial for 



short-range spin glasses in finite dimensions; a competing "droplet" theory (Fisher 



and Huse, 1986 ) approach suggested that many of Parisi's predictions may not sur- 
vive in physical dimensions d = 2,3. The situation is more promising in the case of 
the physically relevant long-ranged Coulomb interaction, where the prediction of the 
mean-field approach may possibly persist even in low dimensions. 

The most interesting test of these ideas relates to the possibility of describing the 
emergence of the "Coulomb gap" in = 2, 3 for localized electrons interacting via 
long-range Coulomb interactions, as first predicted, based on heuristic arguments, a 



long time ago (Efros and Shklovsku, 1975). Here, the single particle density of states 
is predicted to assume a powerlaw form g{e) ~ , notably with a dimensionality- 
dependent exponent 7 = d — 1. 

Conventional mean-field theories typically produce universal dimensionality inde- 
pendent exponents and thus cannot be expected to explain this unfamiliar situation. 



On the other hand, the EDMFT (Chitra and Kotliar, 2000) approach does include the 



effects of spatial correlations, as the bosonic collective modes describing the "cavity 
field" are treated at a Gaussian level, similar to the familiar Hertz-Millis theories of 



quantum criticality (Hertz, 1976 Millis, 1993 ). When applied to the case of long-range 
Coulomb interactions, the corresponding plasmon propagator does reflect (Chitra and 
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Kotliar, 2000 ) both the specific form of the long-range interactions and the dimension- 



ality of the system. When applied to disordered electrons, the replica version of this 
method produces self-consistency conditions which assume a very similar form as in 
ordinary Parisi theory or the simple infinite range electron glass model we examined 
in previous sections. Indeed, recent work has generalized our EDMFT method to the 
long-range models, where the glassy phase of the generalized Coulomb interactions 
(Pankov and Dobrosavljevic, 20051 of the form V{R) ~ has been examined in 

finite dimensions. We will not elaborate of the details of these theories here, but we 
emphasize that most qualitative features of the resulting glassy phase have been found 



to share very similar behavior (Miiller and lofFe, 20041 as for the infinite-range model 



of electron glasses, including the corresponding pattern of replica symmetry break- 



ing (Pankov, 2006 Miiller and Pankov, 2007). Most remarkably, however, this theory 



produced (Pankov and Dobrosavljevic, 2005) an interaction-range and dimensionality 



dependent Coulomb gap exponent 



7 = (d — a)/a, 



(1.95) 



in precise agreement with an appropriate generalization of the Efros-Shklovskii argu- 



ment (Efros and Shklovskii, 1975). This early result was later confirmed (Miiller and 



Pankov, 2007 ) by a more detailed low-temperature solution of the same equations deep 
in the replica-symmetry broken phase, based on a T = new solution of the Parisi 



equation due to Pankov ( Pankov, 2006 ) 



Another striking result of these theories should be noted. Namely, one finds (Pankov 



and Dobrosavljevic, 2005) that the Coulomb pseudogap starts to form around a 
crossover temperatures T* ~ Tc (the Coulomb energy), much above the glass tran- 
sition temperature Tq « 0.05Tc, and in perfect quantitative agreement (Fig. 1.16 



with earlier numerical work (Grannan and Yu, 1993). The EDMFT theory, in fact. 



demonstrated the following: (1) the universal value of the Efros-Shklovskii exponent 
7 is a direct consequence of the marginal stability of the glassy phase; (2) in contrast, 
the emergence of a non-universal Coulomb pseudogap in the high temperature regime 
Tg < T < T* is not related to glassy freezing and the breakdown of screening ( Miiller 



and loffe, 2004 Miiller and Pankov, 2007), and is a robust effect found (Efros, 1992) 



even in the absence of disorder. The two mechanisms for pseudogap formation have 



often been confused ( Grannan and Yu, 1993 ) in previous work, leading to incorrect and 



even misleading interpretations of what determines its form in a given temperature 
range. 

• Quantum melting of the electron glass 

Next, we investigate how the glass transition temperature can be depressed by 
quantum fiuctuations introduced by inter-site electron tunneling. For simplicity, we 
again focus on the simplest infinite-range model of the electron glass. As in other 
quantum glass problems, quantum fluctuations introduce dynamics in the problem, 
and the relevant self-consistency equations cannot be solved in closed form for general 
values of the parameters. In the following, we will see that in the limit of large random- 
ness an exact solution is possible. The main source of difficulty in general quantum 
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Fig. 1.16 The analytical EDMFT predictions ( Pankov and Dobrosavljevic, 2005 I for the sin- 



gle-particle density of states (full lines) are found to be in excellent quantitative agreement 
with simulation results (dashed lines), with no adjustable parameters. Shown are results 



for the three dimensional case studied in Ref. (Sarvestani et al, 19951, corresponding to 
W = l/(2\/3), and temperatures T = 0.4,0.2,0.1,0.05 (top panel), and the two dimensional 
model of Ref. ( |Efros, 1992 1, corresponding to VK = 0, T = 0.1, K = 0.2. Green lines corre- 



spond to the fluid pseudogap phase, while the red curve corresponds to T = 0.05, close to 
the glass transition temperature. 

glass problems relates to the existence of a self-consistently determined "memory ker- 
nel" x(t — r') in the local effective action. By the same reasoning as in the classical 
case, one can also ignore this term since this quantity is also bounded. 

The remaining action is that of noninteracting electrons in the presence of a strong 
random potential. The resulting local compressibility then takes the form 



j3 /■+°° 1 
Xiocis) ^ ^ J dwpe(a;)cosh"^(-/3a;). (1.96) 

Here, Pe(w) is the local density of states, which in the considered large z limit is 
determined by the solution of the CPA equation 
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Fig. 1.17 Phase diagram as a function of quantum hopping t, temperature T and disorder 
strength W. The glass transition temperature To decreases only slowly (as ~ ^/W) in the 
strong disorder limit. In contrast, the critical value of the hoping element to remains finite 
as VK — >■ 00 



Peioj) = IniG(cj); G{uj) = 

TT 



deP{e) 



(1.97) 



In the limit W/t >> 1, it reduces to a narrow resonance of width A = 7r<^P(0) ^ t'^/W 



TT (UJ 



A2 



(1.98) 



The resulting expression for the quantum critical line in the large disorder limit 
takes the form 

taiT = 0,W ^ oo) ^V/y^. (1.99) 

At first glance, this result is surprising, since it means that a, finite value of the Fermi 
energy is required to melt the electron glass at T = 0, even in the W ^ oo limit! This is 
to be contrasted with the behavior of Tq in the classical limit, which according to Eq. 



(1.92 1 was found to decrease as l/W for strong disorder. At first puzzling, the above 



result in fact has a simple physical meaning. Namely , the small resonance width (or 
"hybridization energy") A ~ t'^ /W can be interpreted (Anderson, 1958 Dobrosavljevic 
and Kotliar, 1997) as the characteristic energy scale for the electronic motion. As 
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first pointed out by Anderson (Anderson, 1958), according to Fermi's golden rule. 



the transition rate to a neighboring site is proportional to A and not t, and thus 
becomes extremely small at large disorder. Thus the "size" of quantum fluctuations, 
which replace the thermal fluctuations at T = 0, is proportional to A ~ 1/W, and thus 
becomes very small in t he lar ge W limit. We can now easily underst and t he qualitative 
behavior shown in Eq. ( 1.99 ) by replacing Tq — A ~ ^g/^ ( 1-92 ). The leading 

W dependence cancels out, and we find a finite value for in the W ^ oo limit. 

More generally, we can write an expression for the glass transition critical line in 
the large disorder limit, as a function oi (3 = 1/T and t in the scaling form 



1 = {v/tfmvw), 



with 



iz) - -,z^ 



+ C30 



dx 



+ 00 2^ 

dy 



1 



TT 1 + {x — y)' 



cosh '^{-zy) 



(1.100) 



(1.101) 



At finite disorder an exact solution is not possible, but we can make analytical 
progress motivated by our discussion of the large W limit. Namely, one can imagine 
evaluating the required local compressibilities in Eq. (1.88) by a "weak coupling" 



expansion in powers of the interaction V . To leading order, this means evaluating 
the compressibilities at y = 0, an approximation which becomes exact for W large. 
Such an approximation can be tested for other spin glass problems. We have carried 
out the corresponding computations for the infinite range Ising spin glass model in 
a transverse field, where the exact critical transverse field is known from numerical 
studies. We can expect the leading approximation to underestimate the size of the 
glassy region, i. e. the critical field, since the omitted "memory kernel" introduces long 
range correlations in time, which make the system more "classical" . Indeed, we find 
that the leading approximation underestimates the critical field by only about 30%, 
whereas the next order correction gives an error of less than 5%. Encouraged by these 
arguments, we use this "weak-coupling" approximation for arbitrary disorder strength 
W . Again, the computation of the compressibility reduces to that of noninteracting 
electrons in a CPA formulation; the resulting phase diagram is shown in Fig. |1.17[ 

• Quantum critical behavior of the electron glass 

So far, we have seen how our extended DMFT equations can be simplified for large 
disorder, allowing an exact computation of the phase boundary in this limit. In our 
case, this quantum critical line separates a (non-glassy) Fermi liquid phase from a 
metallic glass phase which, as we will see, features non-Fermi liquid behavior. If one is 
interested in the details of the dynamics of the electrons near the quantum critical line, 
the above simplifications do not apply and one is forced to self-consistently calculate 
the form of the "memory kernel" (local dynamic compressibility) x(t— r'). Fortunately, 
this task can be carried out using methods very similar to those developed for DMFT 



models for metallic spin glasses (Read et ai, 1995). Formulating such a theory is 



technically possible because the exact quantum critical behavior is captured when the 



relevant field theory is examined at the Gaussian level (Miller and Huse, 19931, in the 
considered limit of large dimensions. 
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Because of the technical complexity of this calculation, we only report the main 



results, while the details can be found in (Dalidovich and Dobrosavljevic, 2002). In this 



paper, the full replica-symmetry broken (RSB) solution was found, both around the 
quantum critical line and in the glassy phase. In the Fermi liquid phase, the memory 
kernel was found to take the form 

V'^Xi^n) = D{uJn) + /3(7EA'5u;„,0, 

with 

Here, A is a characteristic energy scale that vanishes on the critical line, which 
also determines the crossover temperature scale separating the Fermi liquid from the 
quantum-critical regime. In contrast to conventional quantum critical phenomena, but 
similarly to metallic spin glasses, the "gap" scale A is equal to zero not only on the 
critical line, but remains zero throughout the entire glassy phase. As a result, the ex- 
citations in this region assume a non-Fermi liquid form 

D{ujn) = -yq^Al^^ - Vl^- 

This behavior reflects the emergence of soft "replicon" modes (iMezard et al., 1986) 



describing in our case low energy charge rearrangements inside the glassy phase. At 
finite temperatures, electrons undergo inelastic scattering from such collective excita- 
tions, leading to a temperature dependence of the resistivity that takes the following 
non-Fermi liquid form 

/9(T) = p(o) + yir3/2. 



Interestingly, very recent experiments (Bogdanovich and Popovic, 2002 ) on two-dimensional 



electron gases in silicon have revealed precisely such a temperature dependence of the 
resistivity. This behavior has been observed in what appears to be an intermediate 
metallic glass phase separating a conventional (Fermi liquid) metal at high carrier 
density from an insulator at the lowest densities. 

Another interesting feature of the predicted quantum critical behavior relates to the 
disorder dependence of the crossover exponent (j> describing how the gap scale A ~ 5r'^ 



vanishes as a function of the distance Sr from the critical line. Calculations (Arrachea 



et al., 2004) show that = 2 in the presence of site energy disorder, which for our 
model plays the role of a random symmetry breaking field, and cj) — 1 in its absence. 
This indicates that site disorder, which is common in disordered electronic systems, 
produces a particularly large quantum critical region, which could be the origin of the 
large dephasing observed in many materials near the metal-insulator transition. 

• Glassy behavior near the Mott- Anderson transition 

As we have seen, the stability of the glassy phase is crucially determined by the 
electronic mobility at T = 0. More precisely, we have shown that the relevant energy 
scale that determines the size of quantum fluctuations introduced by the electrons is 
given by the local "resonance width" A. It is important to recall that precisely this 



quantity may be considered ( Anderson, 1958 ) as an order parameter for Anderson lo- 



calization of noninteracting electrons. Recent work (Dobrosavljevic and Kotliar, 1997 
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Dobrosavljevic and Kotliar, 1998 ) demonstrated that the typical value of this quantity 



plays the same role even at a Mott- Anderson transition. We thus expect A to generally 
vanish in the insulating state. As a result, we expect the stability of the glassy phase to 
be strongly affected by Anderson localization effects, as we will explicitly demonstrate 
in the next section. 



On physical grounds, one expects the quantum fluctuations (Pastor and Dobrosavl 



jevic, 1999) associated with mobile electrons to suppress glassy ordering, but their 



precise effects remain to be elucidated. Note that even the amplitude of such quantum 
fluctuations must be a singular function of the distance to the MIT, since they are 
dynamically determined by processes that control the electronic mobility. 

To clarify the situation, the following basic questions need to be addressed: (1) Does 
the MIT coincide with the onset of glassy behavior? (2) How do different physical 
processes that can localize electrons affect the stability of the glass phase? In the 
following, we provide simple and physically transparent answers to both questions. 
We find that: (a) Glassy behavior generally emerges before the electrons localize; (b) 



while Mott localization (Mott, 1990) tends to suppress it 



Anderson localization (Anderson, 1958) enhances the stability of the glassy phase, 



In order to be able to examine both the effects of Anderson and Mott localization, 
we concentrate on extended Hubbard models given by the Hamiltonian 



Here, Sn^ = rii— {rij) represent local density fluctuations ( (n^) is the site-averaged 
electron density), U is the on-site interaction, and Si are Gaussian distributed random 
site energies of variance . In order to allow for glassy freezing of electrons in the 
charge sector, we introduce weak inter-site density-density interactions Vij, which we 
also choose to be Gaussian distributed random variables of variance /z (z is the co- 
ordination number). We emphasize that, in contrast to previous work (Dobrosavljevic 



et al., 2003a), we shall now keep the coordination number z finite, in order to allow 



for the possibility of Anderson localization. To investigate the emergence of glassy 



ordering, we formally average over disorder by using standard replica methods (Do- 



brosavljevic et al, 20036) and introduce collective Q-fields to decouple the inter-site 
V-teim ( Dobrosavljevic et al., 2003 &| ). A mean-field is then obtained by evaluating the 
Q-fields at the saddle-point level. The resulting stability criterion takes a form similar 
to the one previously discussed (Eq. |1.87[ ) 



i 



0. 



(1.102) 



Here, the non-local static compressibilities are defined (for a fixed realization of disor- 
der) as 

Xij ^ -dni/dsj, (1.103) 

where rii is the local quantum expectation value of the electron density, and [• • ■]dis 
represents the average over disorder. Obviously, the stability of the glass phase is 



determined by the behavior of the fourth-order correlation function x 



(2) 



3 



I dis 
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in the vicinity of the metal-insulator transition. We emphasize that this quantity is 
to be calculated in a disordered Hubbard model with finite range hopping, i.e. in the 
vicinity of the Mott- Anderson transition. The critical behavior of x*^^^ is very difficult 
to calculate in general, but we will see that simple results can be obtained in the limits 
of weak and strong disorder, as follows. 

1.4.3 Large disorder 

As the disorder grows, the system approaches the Anderson transition at t — tc{W) ^ 
W. The first hint of singular behavior of x''^'' in an Anderson insulator is seen by 
examining the deeply insulating, i. e., atomic limit W ^ t, where to leading order we 
set i = and obtain Xij = ^{^i — ^)6ij, i.e. x^'^'^ = [^^^(^i — — +oo! Since we expect 
all quantities to behave in qualitatively the same fashion throughout the insulating 
phase, we anticipate x*^^^ to diverge already at the Anderson transition. Note that, since 
the instability of the glassy phase occurs already at x^^^ = the glass transition 

must precede the localization transition. Thus, for any finite inter-site interaction V ^ 
we predict the emergence of an intermediate metallic glass phase separating the Fermi 
liquid from the Anderson insulator. Assuming that near the transition 

X(^)^:^((VW^)-S)-" (1.104) 



(A and B = tc/W are constants of order unity), from Eq. (1.102) we can estimate the 
form of the glass transition line and we get 

St{W) = tciW) - tc(W) - y2/ap^l-2/a. (1.105) 

The glass transition and the Anderson transition lines are predicted to converge at 
large disorder for a < 2 and diverge for a > 2. Since all the known exponents charac- 
terizing the localization transition seem to grow with dimensionality, we may expect 
a particularly large metallic glass phase in large dimensions. 

Bethe lattices. In order to confirm this scenario by explicit calculations, we com- 
pute the behavior of x*-^^ at the Anderson transition of a half-filled Bethe lattice of 



coordination z = 3. We use an essentially exact numerical approach ( Dobrosavljevic 



and Kotliar, 19971 based on the recursive structure of the Bethe lattice (Abou-Chacra 
et at, 1973). In this approach, local and non-local Green's functions on a Bethe lat- 
tice can be sampled from a large ensemble and the compressibilities Xij can then be 
calculated by examining how a local charge density rii is modified by an infinitesimal 
variation of the local site energy Sj on another site. To do this, we have taken special 
care in evaluating the local charge densities rii by numerically computing the required 
frequency summations over the Matsubara axis, where the numerical difficulties are 
minimized. Using this method, we have calculated x^^-* as a function of W/t (for this 



lattice at half-filling Ep = 2\/2t), and find that it decreases exponentially (Mirlin and 



Fyodorov, 1991 ) as the Anderson transition is approached. We emphasize that only 
a finite enhancement of x^^^ is required to trigger the instability to glassy ordering, 
which therefore occurs well before the Anderson transition is reached. The resulting 
r = phase diagram, valid in the limit of large disorder, is presented in Fig. |1.18[ 
Note that the glass transition line in this case has the form tQ{W) ~ W, in agreement 
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Fig. 1.18 Phase diagram for tlie z = 3 Betlie lattice, valid in the large disorder limit. The 
inset shows x'^^ as a function of disorder W. 

with the fact that exponential critical behavior of x*-^-* corresponds to a — >■ cx) in the 
above general scenario. These results are strikingly different from those obtained in a 



theory which ignores localization (Pastor and Dobrosavljevic, 1999), where tciW) was 
found to be weakly dependent on disorder and remain finite as W — > oo. Anderson 
localization effects thus strongly enhance the stability of the glass phase at sufficiently 
large disorder. Nevertheless, since the Fermi liquid to metallic glass (FMG) transition 
occurs at a finite distance before the localization transition, we do not expect the 
leading quantum critical behavior ( Dalidovich and Dobrosavljevic, 2002[ ) at the FMG 
transition to be qualitatively modified by localization effects. 

Typical medium treatment. As an alternative approach to the Bethe lattice calcu- 
lation, in this section we introduce Anderson localization to the problem by using the 
formaHsm of "Typical Medium Theory" (TMT) ([Dobrosavljevic al, 2003a[), which 



was explained in detail in Section 1.2 We calculate the cavity field AxYpit^) by solv- 
ing the relevant self-consistency condition (Dobrosavljevic et at, 2003a), which in turn 
allows us to find local compressibilities: 



dn 
de. 



l_d_ 
TT de^ 
1 



dujImG{ei, uj, W), 



- Si - A 



TYP 



(1.106) 
(1.107) 



needed to determine the critical line of the glass transition. These calculations were 
performed using a model semicircular bare density of states po(w) and a box dis- 



tribution of disorder P{si). The resulting phase diagram is shown in Fig. 1.19 The 
intermediate metallic glassy phase still exists, but shrinks as — >■ oo, reflecting the 
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Fig. 1.19 Phase diagram from Typical Medium Theory of Anderson locahzation (Dobrosavl- 
jevic et al, 2003ttl, giving a = 1. The intermediate metalhc glass phase shrinks as disorder 



W grows. Compare this to the Bethe lattice case (Fig. 1.181, where a = oo 



small value of the critical exponent a — 1, vifhich can be shown analytically within 
TMT. A more realistic value for this exponent, corresponding to d = 3, requires more 
detailed numerical calculations, which remains a challenge for future work. 

Low disorder - Mott transition. In the limit of weak disorder W -^U^V ^ interactions 
drive the metal-insulator transition. Concentrating on the model at half-filling, the 



system will undergo a Mott transition (Mott, 1990) as the hopping t is sufficiently 
reduced. Since for the Mott transition tMott{U) ^ U, near the transition W t and 
to leading order we can ignore localization effects. In addition, we assume that V <^ U 
and to leading order the compressibilities have to be calculated with respect to the 
action Sei of a disordered Hubbard model. The simplest formulation that can describe 
the effects of weak disorder on such a Mott transition is obtained from the dynamical 
mean-field theory (DMFT) (Georges et al, 19961. This formulation, which ignores 
localization effects, is obtained by rescaling the hopping elements t — > t/y/z and then 
formally taking the limit of large coordination z ^ oo (see Section [1.1.4 ) . To obtain 



qualitatively correct analytical results describing the vicinity of the disordered Mott 
transition at T = 0, we have solved the DMFT equations using a 4-boson method 



(Dobrosavljevic et al., 20036). At weak disorder, these equations can be easily solved 
in closed form and we simply report the relevant results. The critical value of hopping 
for the Mott transition is found to decrease with disorder as 



t,{W)^t°Jl~A{W/U)^ + ■■■), 



(1.108) 



where for a simple semi-circular density of states (Georges et al., 1996) t° = SttU/QA 
(in this model, the bandwidth B — it). Physically, the disorder tends to suppress the 
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Mott insulating state, since it broadens the Hubbard bands and narrows the Mott- 
Hubbard gap. At sufficiently strong disorder W >U, the Mott insulator is suppressed 
even in the atomic limit t — > 0. The behavior of the compressibilities can also be 
calculated near the Mott transition and to leading order we find 



X 



(2) _ 



(1- 



t ' 



{l + 2S{W/Uf). 



(1.109) 



Therefore, as with any compressibility, x'^^ is found to be very small in the vicinity of 
the Mott transition, even in the presence of finite disorder. As a result, the tendency 
to glassy ordering is strongly suppressed at weak disorder when one approaches the 
Mott insulating state. 

Finally, having analyzed the limits of weak and strong disorder, we briefly comment 
on what may be expected in the intermediate region W ~ ?7. The Mott gap cannot 
exist for W > U , so in this region and for sufficiently small t (i. e. kinetic energy), 
one enters a gapless (compressible) Mott- Anderson insulator. For W ^ U, the compu- 
tation of x^^^ requires the full solution of the Mott- Anderson problem. The required 



calculations can and should be performed using the formulation of (Dobrosavljevic 



and Kotliar, 1997), but that difficult task is a challenge for the future. However, based 
on general arguments presented above, we expect x^"^^ to vanish as one approaches the 
Mott insulator {W < U), but to diverge as one approaches the Mott- Anderson insu- 
lator {W > U). Near the tetracritical point, we may expect x''^'' ~ 5W~^^5t°'^ where 
5W = W — WMott{t) is the distance to the Mott transition line, and 5t = t — tc{W) 



is the distance to the Mott- Anderson line. Using this Ansatz and Eq. (1.102), we find 
the glass transition line to take the form 



St = taiW) - t^{W) - dW^/"; W ^W, 



Mott- 



(1.110) 



We thus expect the intermediate metallic glass phase to be suppressed as the disor- 
der is reduced, and one approaches the Mott insulating state. Physically, glassy behav- 
ior of electrons corresponds to many low-lying rearrangements of the charge density; 
such rearrangements are energetically unfavorable close to the (incompressible) Mott 
insulator, since the on-site repulsion U opposes charge fluctuations. Interestingly, very 
recent experiments on low density electrons in silicon MOSFETs have revealed the 
existence of exactly such an intermediate metallic glass phase in low mobility (highly 
disordered) samples (Bogdanovich and Popovic, 2002). In contrast, in high mobility 
(low disorder) samples (Jaroszyhski et ai, 2002), no intermediate metallic glass phase 
is seen and glassy behavior emerges only as one enters the insulator, consistent with 
our theory. Similar conclusions have also been reported in studies of highly disordered 



In02 films (Ben-Chorin et al., 1993 Ovadyahu and Pollak, 1997), where the glassy 



slowing down of the electron dynamics seems to be suppressed as the disorder is re- 
duced and one crosses over from an Anderson-like to a Mott-like insulator. In addition, 



these experiments (Bogdanovich and Popovic, 2002 Jaroszyhski et al., 2002) provide 



striking evidence of scale-invariant dynamical correlations inside the glass phase, con- 
sistent with the hierarchical picture of glassy dynamics, that generally emerges from 
mean-field approaches (Mezard et al, 1986) such as the one used in this work. 
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1.5 Beyond DMFT: loop expansion and diffusion modes 

DMFT approaches to correlation and disorder have already provided significant insight 
in a number of key phenomena and processes around the metal-insulator transition. 
Most importantly, these methods have made it possible to describe strong correla- 
tion effects associated with both "Mottness" and the glassy behavior of electrons. 
The information provided is, in many ways, rather complementary to the conventional 
weak-coupling approaches which focused on the effects of long-wavelength diffusion 
modes within the Fermi liquid picture. Because of their focus on local correlation ef- 
fects, the DMFT theories are ill-suited to describe those phenomena and regimes where 
long-range spatial correlations dominate. One of the biggest challenges for future work 
is to find ways to combine the strength of both approaches and introduce systematic 
methods to incorporate the non-local spatial correlations ignored by DMFT. While 
little has been done so far to implement this ambitious program, some approaches 
have already been outlined that provide guidance on how the problem should be for- 
mulated. We conclude our discussion by a brief outline of a promising approach to 
include both the local correlation processes and the diffusion-mode effects within a 
single theory. 

1.5.1 Gauge Invariant Models of Wegner 

The task of identifying the leading nonlocal corrections to DMFT can be formulated 
in a particularly elegant and transparent fashion for certain special models of disor- 



der. Here, the DMFT equations can be obtained ( Dobrosavljevic and Kotliar, 1994) 
from a functional-integral formulation at the saddle-point level, allowing systematic 
corrections using a loop expansion. In this class, the inter-site hopping elements are 
assumed to be random variables of the form 

tij = Vij 9{xi,Xj), (1-111) 
and in addition, there can be an arbitrary distribution of site energies Ej. Here, Uij-s 
are independent bond variables with a symmetric distribution, i. e. vlj''^^ = 0, and 
g{xi,Xj) is an arbitrary function of local site variables Xi. 

The special class of models which have a symmetric distribution of hopping ele- 



ments has a very simple physical interpretation. As first observed by Wegner (Wegner, 



1979 Schaffer and Wegner, 1980), in these "gauge invariant" models, the phases of 
the electrons undergo random shifts at every lattice hop and so the mean free path £ 
reduces to one lattice spacing. On general grounds, on length scales longer than £ the 
details of the lattice structure are washed out by disorder, so that for gauge invariant 
models in the large dimensionality limit, the details of the lattice structure become 
irrelevant. We contrast this with the models with arbitrary disorder discussed previ- 
ously, which have a well defined pure limit and accordingly can also have an arbitrarily 
large mean free path. The presence of this intermediate length scale {£ can be much 
larger than the lattice spacing a, but much smaller than the localization length ^) 
is often irrelevant to both long-wavelength phenomena such as localization and local 
phenomena such as the Mott transition. The gauge invariant models avoid these un- 
necessary complications without disrupting any of the qualitative properties on either 
very short or very long length scales. 
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The general properties are the same for all the models in this class, but for the 
simplicity of our presentation, we will restrict our attention to the separable case 
( Shiba, 1971 1 where g{xi, xj) = XiXj , with an arbitrary distribution Px {x. 
variables Xi. For the trivial choice of Px{xi) 



for the site 
1), the models reduce to the gauge 



invariant models of Wegner ( Wegner, 1979[ Schaffer and Wegner, 1980). Nontrivial 
distributions Px{xi) which extend to small values of the variable Xi are useful for 
the study of disorder-induced local moment formation (Milovanovic et al., 1989). Sites 
with small Xi represent sites with weak hybridization. At intermediate correlations, we 
expect sites with small Xi to behave as local moments and give large contributions to 
thermodynamic quantities such as the specific heat coefficient 7 = C/T, while other 
sites remain in the itinerant regime. 

Also, we take the ytj-s to be Gaussian random variables with zero mean and vari- 
ance 

1 



2 _ 



- fij t 



Here, the (uniform) matrix fij specifies the lattice structure 



1, for connected sites, 
0. for disconnected sites. 



(1.112) 



(1.113) 



and we have scaled the (square of the) hopping elements by the coordination number 



in order to obtain finite results in the z — ?► 00 limit. 



1.5.2 Functional Integral Formulation 

At this point, it is convenient to explicitly perform the averaging over the Gaussian 



random (bond) variables y^j, using the standard replica formulation (Wegner, 1979 
[Schaffer and Wegner, 1980|). The hopping part of the action then assumes the form 
(Dobrosavljevic and Kotliar, 1994) 



Si 



hop 



fij Xi X 



dr[clt{T)cZJr) + h.c: 



(1.114) 



As we can see from this expression, the averaging over disorder has generated a quartic 
term in the action that is nonlocal in (imaginary) time, spin and replica indices. We 
are now in a position to introduce collective Q-fields (Wegner, 1979 Schaffer and 
Wegner, 1980{|Finkerstein, 1983{|Finker stein, 1984) of the form (in terms of Matsubara 
frequencies ui = 2mTT; the indices "n" are omitted for brevity) 



1 



E 



2 t " 



(1.115) 



by decoupling the (quartic) hopping term using a Hubbard-Stratonovich transforma- 
tion. For simplicity, as before, we will ignore the superconducting phases, as well as 
the fluctuations in the particle-particle (Cooper) channel, so that the Q-field does not 
have anomalous components. The procedure can be straightforwardly generalized to 
include the omitted terms (Efetov et al, 1980). 
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It is now possible to formally integrate out the electron (Grassmann) fields and the 
resulting action for the Q-helds can be written as 

S[Q] = ShopiQ] + SiociQ]- (1.116) 

The nonlocal part of the action ShopiQ] takes a simple quadratic form in terms of the 
Q fields 

ShoAQ] = E E E E QZlZr''H^)QZlZl''''Hj), (1-117) 

ij aia2 siS2 (^1(^2 

where, Kij = \f^^ is the inverse lattice matrix, scaled by coordination number z. In 
contrast, all the nonlinearities are contained in the local part of the action 

S'/oc[Q] = -E j dxiPxixi) j dSiPs^Si) j Dc\Dcieyi■p^~Seff[c\,Ci,Q^,x^,ei]^, 

(1.118) 

where the local effective action takes the form 

- E E E 4"/(^i)[(*^i +M-e.)5o,o.<5.,..<s.,.. -X? t'Qz\zv'^''{^]cZM) 

OL\OL2 S1S2 UJ1UJ2 

+^E E ci;(^i)<t(^2)4;(a.3)c^;(^4). (i.iw) 

The local effective action Seff[cl,Ci,Qi,Xi,ei] is identical to the action of a (gener- 
alized) Anderson impurity model embedded in an electronic bath characterized by a 
hybridization function t'^ QZ]ZV''"' ^ i'i-)- We can thus interpret our system as a col- 
lection of Anderson impurity models (Anderson, 1961) that are "connected" through 
the existence of collective Q-fields. Here we note that, in contrast to an ordinary An- 
derson model, the hybridization function is now non-diagonal in frequency, spin and 
replica indices. Physically, this reflects the fact that in general dimensions a given site 
can be regarded as an Anderson impurity model in a fluctuating bath which breaks 
translational invariance in time, space and spin. 

1.5.3 Saddle-Point Solution 

The action has a general form which is very similar to standard lattice models inves- 
tigated in statistical mechanics (Goldenfeld, 1992). As usual, the problem simplifies 
considerably in the limit of large coordination number, when the spatial fluctuations 
of the Hubbard-Stratonovich field {Q in our case) are suppressed and the mean-field 
theory becomes exact. It is worth pointing out that there are two classes of lattices 
which can have large coordination: 

(a) Lattices with short-range bonds but living in a space of large dimensionality. 
For example, on a hypercubic lattice with nearest neighbor hopping in d dimensions, 
z = 2d. 
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(b) Lattices embedded in a finite dimensional space but having long hopping range. 
In this case, the lattice matrix fij takes the form 

■^'^ ^ { 0, otherwise, (1-120) 

and the coordination number z ^ L"^. 

In either case, when z — > oo, the functional integral over Q-fields, representing the 
partition function, can be evaluated (exactly) by a saddle-point method and we obtain 
a mean-field theory. In order to derive the mean-field equations in our case, we look 
for extrema of the action S[Q] with respect to variations of the Q-fields, i. e. 

^^[Q] -0 (1121) 

Since the saddle-point solution is translationally invariant in time and space and con- 
serves spin, it is diagonal in all indices 

[QZlZ'/'^Hi)] \sp = S^,^Js,sJu:,u.,Q's''H, (1-122) 
and the saddle-point equations assume the form 



(cj) = / de^Psie^) / dx,Px{x,)x^^ G,,,(w), (1.123) 



where 

If we identify 



G^.s{uJ) =< cl{uj)Cs{u!) >S,„ [ct x.QS^- ,Xi ■ (1.124) 



WsA.uj) = x^,t^Q^/{Lj), (1.125) 



we see that our saddle-point equations become identical to the standard DMFT (d — > 
oo) equations, when applied to the appropriate model of hopping disorder. We empha- 
size that the present equations are exact at z ^ oo for an arbitrary lattice, due to the 
presence of the "gauge invariant" form of the hopping disorder. Since the saddle point 
equations determine the local effective action, this means that all the local correlation 
functions will be insensitive to the lattice structure in this mean-field limit. However, 
other properties such as the tendency to the formation of spin and charge density 
waves are very sensitive to the details of the lattice structure. 

As an example, we can compare the case of simple hopping disorder, t.^ — yij, in a 
bipartite lattice such as the Bethe lattice and the case of a lattice with infinite range 
hopping (the limit L — cx) of the model (b) above). The self-consistency (mean-field) 
equations are identical in the two cases and in fact reduce to those of a pure Hubbard 
model on a Bethe lattice with hopping t. On the other hand, it is well established 



(Rozenberg et al., 1992) that in the first case the system is unstable towards the 
formation of an antiferromagnetic ground state, even for arbitrarily small U/t, while in 
the second, the system remains paramagnetic for any U/t, due to the large frustration. 



In many physical systems, such as doped semiconductors (Paalanen and Bhatt, 1991 ), 
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disorder introduces large amounts of frustration and magnetic ordering does not occur, 
even though the system is strongly correlated. In order to study such situations, it is 
useful to have at hand microscopic models that have a non-magnetic ground state 
and allow one to study the approach to the metal-insulator transition which occurs at 

r = o. 

1.5.4 Loop Expansion 

The present approach is particularly convenient for the study of the effects of strong 
correlations on disorder-driven transitions and the interplay of Anderson localization 
and strong correlations in general. This is especially true since Anderson localization is 
not present in d = cx) (or infinite range) models and so one has to extend the approach 
to include spatial fluctuations missing from the mean-field description. In order to 
systematically study the fluctuation effects, we proceed to carry out an expansion in 
terms of the deviations of the collective Q-fields from their saddle-point value, i. e., 
in powers of SQ{i) — Q{i) — Q^^ . This procedure, also known as a loop expansion 



(Goldenfeld, 1992) has been used in other disordered problems, such as spin-glasses 
( de Dominicis, C. et al., 1991[ ), to generate systematic corrections to the mean-field 



theory. The method is particularly convenient when applied to long-range models (de 



Dominicis, C. et al., 1991) (class (b) above), since in that case the loop corrections 
are ordered by a small parameter 1/z. The loop expansion can be applied also to 
large dimensionality models, but in that case a given order in a loop expansion can be 



considered to be an infinite resummation of the simple 1/d expansion (Georges et al, 



1990), since each term contains all powers of 1/d. 

When the expansion of the effective action in terms of 5Q-s is carried to lowest, 
quadratic order, we obtain a theory describing gaussian fluctuations around the sad- 



dle point, which represent weakly interacting collective modes (Wegner, 1979 Schaffer 



and Wegner, 1980 Finkel'stein, 1983 Finkel'stein, 1984). Higher order terms in the 



expansion then generate effective interactions of these modes, which under appropri- 
ate conditions can lead to fluctuation-driven phase transitions. In practice, if all the 
components of the collective Q-fields are retained in this analysis, the theory becomes 
prohibitively complicated and cumbersome. However, in order to analyze the critical 
behavior, it is not necessary to keep track of all the degrees of freedom, but it suffices 
to limit the analysis to soft modes, i. e., those that represent low energy excitations. 
In disordered metallic phases, charge and spin conservation laws lead to the existence 
of diffusion modes, which are the hydrodynamic modes describing charge- and spin- 



density relaxation. In the Fermi liquid regime (Castellani et al., 1987), all the other 
collective excitations require higher energy and can be neglected in a hydrodynamic 
description of the system. One is then led to construct a theory of interacting diffusion 
modes, as a theory of critical phenomena for disordered interacting systems. 

This line of reasoning was used in field-theoretical approaches to the localization 



problem of noninteracting electrons, as first developed by Wegner (Wegner, 1979 



Schaffer and Wegner, 19801. In this theory, collective Q-fields, similar to the ones 
presented in this paper, are introduced. At the saddle-point level, no phase transition 
occurs and all the states are extended. An analysis of the fluctuations of the Q-Helds 
is then performed and a subset of those fluctuations identified, which represent the 
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hydrodynamic (diffusion) modes. Only the fluctuations of these modes are retained and 



an effective hydrodynamic theory is constructed - the non-linear sigma model ( Wegner, 



1979 Schaffer and Wegner, 1980 1 . The interactions of these modes lead to the metal- 



insulator (localization) transition, which was analyzed using renormalization-group 



techniques and 2-\-e expansions. In subsequent work, Finkelshtein ( Finkel'stein, 1983 



Finkel'stein, 1984) was able to apply a similar procedure to interacting disordered 
electrons. However, this theory is based on a number of implicit assumptions that 
restrict its validity to Fermi-liquid regimes. In the language of Q-fields, this theory 
again expands around a noninteracting saddle-point and the interaction effects appear 
only at the level of the fluctuation corrections. 

When strongly correlated electronic systems are considered, much of the physics 
relates precisely to the destruction of coherent quasiparticles by inelastic processes, 
so that one needs a description that is not limited to Fermi-liquid regimes. In the 
language of hydrodynamics, new soft modes appear, which indicate the tendency to- 
wards interaction-driven instabilities. In particular, strong correlations can lead to 
local moment formation and the Mott transition. In both cases, charge fluctuations 
are suppressed and low-energy spin fluctuations dominate the physics. 

In the present approach, in contrast to the work of Wegner and Finkelshtein, the 
strong correlations are treated in a non-perturbative fashion already at the saddle- 
point (mean-field) level, so that all the soft modes can be systematically included. 
In particular, even at the one-loop level, we can address the question of how the 
disorder-induced local moment formation and the approach to the Mott transition 
affect the weak localization (diffusion) corrections. Ultimately, our approach indicates 
how a more general low energy theory can be constructed that extends the (T-model 
description so as to include strong correlation effects. 

In the present discussion, we will limit our attention to the form of the Gaus- 
sian fluctuations of the Q-flelds, which allow one in principle to compute the leading 
corrections to mean-field theory. The Gaussian part of the action takes the form 

-t'W{W)W{h)5i,i,5i,u + t^T{h ■ ■ ■ k)] SQi,u{-k). (1.126) 

This expression is appropriate for the long-ranged model (b) above, in which case the 
inverse lattice matrix in momentum space takes the form K{k) « 1 -I- L'^k^ and we 
cut off the momentum integrals at A = 2t:/L. Note that the coefficient of A:^, which 
can be interpreted as the stiffness of the 5Q modes, is ~ L^, so we see that indeed the 
fluctuations are suppressed at L — > oo. In the above formula, the index Im is used to 
represent the frequency, spin and replica indices. The local vertex function r(Zi • • • ^4) 
is given by 



T{li - ■ -li) = j deiPsisi) J dxtPx{xi)xf < cl{li)ct{l2)cl{l3)ci{U) >s^„[Qsr] . 

(1.127) 

At this level, the dynamics of the collective fluctuations 6Q is governed by the form of 
5^^-'[5(5]j which is expressed in terms of the local correlation functions of the saddle- 
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point theory, i. e., of the d = oo disordered Hubbard model. Accordingly, a detailed 
study of the d = oo limit does not provide only a mean-field description of the problem, 
but also determines the form of the leading corrections resulting from fluctuations. 
Extensions of the theory to include the effects of these fluctuation corrections remain 
to be addressed in more detail in future work. 
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